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ABSTRACT 

A  computer-aided  data  acquisition  system  was  developed 
and  a  microthermocouple  probe  constructed  to  obtain  thermal 
distributions  in  turbulent  buoyant  jets  exposed  to  a  cross¬ 
flowing  ambient  fluid.  The  system  performed  high  speed 
temperature  measurements  as  a  microthermocouple  probe  was 
automatically  traversed  through  a  sequence  of  preprogrammed 
positions  under  the  control  of  a  microcomputer.  Operability 
of  the  apparatus  was  demonstrated  by  measuring  temperature 
distributions  in  planes  perpendicular  to  the  streamwise  axis 
of  jets  from  which  contour  plots  of  temperature  were  gener¬ 
ated.  Using  temperature  distributions  along  with  velocity 
distributions  allow  buoyant  jet  characteristics  to  be 
computed,  including  the  entrainment  rate  of  ambient  fluid, 
jet  trajectory,  and  heat  transfer  to  the  ambient.  The 
experimental  technique  is  discussed  and  temperature  contour 
plots  for  a  jet  at  various  planes  are  presented. 
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NOMENCLATURE 


Incremental  Cross-Sectional  Area  in  the 
Temperature  Matrix 

Jet  Half-width 

Normalized  Jet  Half-width 

Specific  Heat 

Diameter  of  the  Jet  at  the  Nozzle 
Binary  Mass  Diffusion  Coefficient 
Densiometric  Froude  Number 
Acceleration  of  Gravity 

Heat  Transfer  Rate  from  the  Jet  to  the  Ambient 
Fluid 

Ambient- to-Nozzle  Flow  Ratio 
Length  of  the  Probe  Arm 

Radial  Distance  from  the  Center  of  the  Jet 
Length  of  the  Probe 
Schmidt  Number 

Streetwise  Coordinate  Along  the  Jet  Centerline 
Normalized  Jet  Temperature 
Ambient  Fluid  Temperature 

Jet  Temperature  Within  the  Temperature  Matrix 
Centerline  Jet  Temperature 
Nozzle  Temperature 

Jet  Temperature  As  Measured  by  the  Probe 

Temperature  Within  the  Jet  at  a  Radial  Distance 
(r)  from  its  Center 


Centerline  Velocity  of  the  Jet  at  the  Nozzle 
Discharge  Velocity  of  the  Jet 

'  .  locity  within  the  Jet  at  a  Radial  Distance  "r" 
t rom  its  Center 

Jet  Velocities  Corresponding  to  Locations  within 
the  Temperature  Matrix 

Normalized  Centerline  Velocity 

Entrainment  Coefficient;  Offset  Angle  of  the 
Probe  Arm 

Offset  Angle  of  the  Probe  Mounting  Bracket 

Probe  Angle  of  Deflection  from  Horizontal 

Local  Angle  of  Inclination  from  Horizontal  of  the 
Jet  Streamwise  Axis 

Spreading  Ratio 

Kinematic  Viscosity 

Density  of  the  Jet  Fluid 

Density  of  the  Ambient  Fluid 

Density  of  the  Ambient  Fluid  at  the  Nozzle  Exit 

Density  of  the  Jet  at  Centerline 

Centerline  Density  of  the  Jet  at  the  Nozzle 

Angle  of  Inclination  of  the  Data  Plane  from  Horizontal 


I .  INTRODUCTION 

Buoyant  jets  are  very  common  in  nature.  We  see  them  in 
the  form  of  exhaust  gases  emitted  from  smoke  stacks  of 
refineries,  mills  and  ships.  We  see  them  in  the  form  of 
heated  waste  water  expelled  into  the  sea  from  power  plants 
and  from  the  main  propulsion  condensers  in  steam  driven 
ships  and  submarines.  It  is  no  wonder  that  the  fluid 
mechanics  and  heat  transfer  characteristics  of  buoyant  jets 
have  been  of  interest  to  environmental,  civil  and  mechanical 
engineers  for  decades.  To  evaluate  their  ecological  impact, 
and  of  most  recent  interest,  to  harness  buoyant  jets  as  a 
means  of  detecting  military  targets  and  guiding  weapons,  it 
is  necessary  to  develop  models  which  accurately  predict 
their  trajectory  and  decay. 

Most  studies  to  date  have  dealt  with  buoyant  jets  rising 
through  a  quiescent  ambient  fluid;  however,  in  nature  most 
problems  involve  flowing  ambient  fluids.  Relatively  little 
experimental  work  has  been  done  with  buoyant  jets  in  cross- 
flow  and,  according  to  Hilder  [Ref.  1],  the  trajectories  of 
jets  and  the  entrainment  rates  of  ambient  fluid  predicted 
from  previous  work  do  not  agree  well  with  one  another.  Most 
mathematical  models  of  buoyant  jets  in  crossflow  assume 
Gauss i an  profiles  for  velocity  and  temperature.  Nickodem 
[ Ref .  2 1  has  shown  through  experiments  that  in  fact,  the 


Gaussian  profiles  of  velocity  are  altered  by  crossflow. 
This  leads  one  to  suspect  that  the  same  may  be  true  for  th 
temperature  profiles. 

The  objective  of  this  work  was  to  develop  a  system  to 
thermally  map  a  buoyant  jet  in  crossflow.  Then,  by  measur 
ing  both  velocity  and  temperature  distributions,  improved 
computations  of  entrainment,  trajectory  and  heat  transfer 
characteristics  of  jets  can  be  made  thereby  giving  rise  to 
more  accurate  models. 


BUOYANT  JETS  DISCHARGED  TO  A  CROSSFLOW 


A.  PROPERTIES  OF  BUOYANT  JETS 

A  buoyant  jet  is  characterized  by  a  momentum  and  a 
density  differential  between  the  jet  and  its  surrounding 
ambient  resulting  from  a  variation  in  temperature  and/or 
fluid  concentrations .  Therefore,  fluid  motion  in  the  jet 
is  governed  by  both  inertial  and  buoyant  forces.  The  non- 
dimensional  ratio  of  these  forces,  known  as  the  densiometric 
Froude  number,  provides  an  important  quantitative  measure¬ 
ment  of  jet  characteristics  and  is  shown  below. 


U 


F  = 


9D(pa  _po)/0o 


where  Uq  is  the  jet's  discharge  velocity,  g  is  the  accelera¬ 
tion  of  gravity,  D  is  the  discharge  diameter  of  the  jet, 
is  the  density  of  the  crossflowing  ambient  and  o  is 
the  density  of  the  jet  fluid  at  its  point  of  discharge. 

The  Gaussian  velocity  and  temperature  profiles  assumed 
by  most  models  of  buoyant  jets  are  very  similar.  Velocity 
behavior  is  given  by: 

U(r)  =  U  exp(-r^/B^) 
m  r 


where 
va  r  i  ab 


U 

m 

le 


is  the  centerline  velocity,  r  is  the  independent 
and  a  radial  distance  from  the  centerline  of  the 
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A  s  r 


jet,  and  B  is  defined  as  nominal  ;ot  ha  1  fw  idth . 
approaches  B,  velocity  decays  to  {  1 , 'e  '  [  Rot  .  3  i  .  Simi¬ 

larly,  temperature  behavior  is  giver,  by: 

T  (  r )  =  T  exp  (  -  r  2  -2B2) 

m 

where  Tm  is  the  centerline  temperature,  r  and  B  are  defined 
the  same  as  above  and  1,  a  spreading  ratio,  is  the  inverse 
of  the  turbulent  Schmidt  number  (S).  S  is  defined  as  the 
ratio  of  the  molecular  momentum  and  mass  dif fusivities  and 
is  equal  to  v/D  where  v  is  the  kinematic  viscosity  and 
D  is  the  binary  mass  diffusion  coefficient  associated  with 
substances  A  and  B  [Ref.  4].  Although  A  varies  inversely 
with  the  Froude  number,  the  change  is  very  slight,  and  in 
the  case  where  substances  A  and  B  are  both  water,  A  is 
slightly  greater  than  1.  Hirst  [Ref.  3]  found  A  to  vary 
between  1.16  at  F  =  0  to  1.11  at  F  -  infinity.  The  net 
effect  then,  is  a  more  gradual  temperature  decay  than  was 
found  with  velocity. 

Most  buoyant  jet  models  consider  the  entrainment  of  the 
ambient  fluid  into  the  jet  and  are  based  on  relevant  conser¬ 
vation  equations  of  mass,  momentum  and  energy.  In  conserva¬ 
tion  of  mass,  the  downstream  change  in  total  mass  of  the  jet 
is  equated  to  the  mass  of  the  entrained  fluid.  The  conser¬ 
vation  of  momentum  must  consider  both  vertical  and  horizontal 
contributions.  Changes  in  vertical  momentum  are  equated 
to  the  buoyant  forces  while  changes  in  the  horizontal 


momentum  of  the  jet  are  equated  to  the  horizontal  momentum 
of  the  entrained  fluid.  The  conservation  of  energy  involves 
energy  changes  resulting  from  variations  in  the  ambient 
temperature  as  caused  by  the  jet.  Hilder  [Ref.  1]  developed 
the  following  governing  equations  in  non-dimensional  differ¬ 
ential  form. 


CONTINUITY 


HORIZONTAL  MOMENTUM 


VERTICAL  MOMENTUM 


ENERGY 


d  ,  .2. 

-5—  u  b 
ds  m 


2  tb  [  u  -R  cos  -  +  a,Rsin 

m  3 


(u2b2  cos  3) 
ds  m 


4Rab[]u  -Rcosr  +a1Rsinv] 
m  1  3 


d  ,  2,2  . 

v—  (u  b  sir.  ?) 
ds  m 


:  -0  ,  .  2.  2 

a  m  ,  2.n  b 


-) 


ao  "mo  F 


2  2 

-  (u_T._Lb  , 


Lt 


ds  m  {,2+1) 


a  -  ,2 

•  u  b  ; 


ds  m 


B.  FLOW  REGIMES 

The  jet  passes  through  several  regimes  as  it  travels 
from  the  nozzle  through  the  ambient.  The  three  regions  most 
frequently  referred  to  are  shown  in  Figure  1.  They  are  the 
zone  of  flow  establishment,  the  zone  of  established  flow 
and  the  far-field  zone  [Ref.  3].  In  the  zone  of  flow  estab¬ 
lishment,  the  velocity  and  turbulence  profiles  transform  from 
the  conditions  within  the  nozzle  to  a  free  turbulent  flow 
condition.  It  is  in  this  region  that  the  jet  begins  to  mix 


with  the  ambient  fluid;  however,  the  flow  is  still  more 


Lgure  1.  Typical  Flow  Regions  in 


strongly  influenced  by  the  nozzle  discharge  conditions  than 
by  the  ambient.  When  the  turbulent  mixing  has  reached  the 
centerline  of  the  jet,  the  zone  of  established  flow  is  said 
to  begin.  In  this  region,  the  profiles  have  assumed  their 
free  turbulent  shapes.  Now  the  flow  is  governed  by  the 
jets'  momentum  and  buoyancy  as  well  as  by  the  condition  of 
the  crossflow.  The  far  field  zone  is  defined  as  that  region 
in  which  jet  momentum  is  depleted  and  the  jet  fluid  is  con- 
vected  and  diffused  by  the  ambient  currents  and  turbulence. 

C.  EFFECTS  OF  CROSSFLOW 

At  the  immediate  exit  of  a  cylindrical  nozzle,  a  verti¬ 
cally  discharged  buoyant  jet  has  a  nearly  uniform  velocity 
distribution  and  has  the  same  cross-sectiona 1  shape  as  the 
nozzle  itself.  The  velocity  gradient  between  the  jet  and 
the  crossflowing  ambient  creates  longitudinal  shear  stresses 
at  the  jet's  sides,  a  positive  pressure  region  immediately 
upstream  and  a  negative  pressure  region  immediately  downstream 
of  the  jet.  This  results  in  the  deflection  of  the  jet's 
trajectory  in  the  downstream  direction  (Figure  2),  the 
creation  of  counterrotating  vortices  at  the  jet's  outer 
edges  and  the  deformation  of  the  original  circular  cross- 
sectional  shape  into  the  form  of  a  kidney.  As  the  stream- 
wise  axis  of  the  jet  approaches  the  direction  of  the  cross- 
flow,  these  effects  become  progressively  less  pronounced. 


17 


Figure  2.  Coordinate  System  for  a  Buoyant  Jet 


Ill . 


EXPERIMENTAL  APPARATUS 


A.  SYSTEM  OVERVIEW 

A  surplus  milling  machine  was  configured  with  synchron¬ 
ous  drive  motors  interfaced  with  a  microcomputer  that  auto¬ 
matically  positioned  its  bed.  It  was  used  as  a  three-dimensional 
positioning  platform  in  the  same  manner  as  in  the  laser 
Doppler  velocimetry  work  undertaken  by  Nickodem  [Ref.  2] . 

The  milling  machine  was  placed  adjacent  to  a  rectangular 
plexiglass  flume  through  which  the  crossflowing  ambient 
fluid  flowed.  A  vertical  nozzle  was  installed  in  the  base 
of  the  flume  to  provide  the  jet.  A  temperature  probe  was 
suspended  through  an  opening  in  the  top  of  the  flume  above 
the  nozzle  by  an  arm  attached  to  a  base  mounted  on  the 
milling  machine  bed  as  shown  in  Figure  3.  As  the  probe  was 
automatically  traversed  through  a  series  of  preprogrammed 
positions  across  the  jet,  temperature  data  was  automatically 
sensed  and  stored  at  high  speeds  by  the  computer. 


B.  CROSSFLOW  SYSTEM 

As  illustrated  in  Figure  4,  the  crossflow  circulation 
pump  took  water  from  the  cylindrical  248.8  1  (65.7  gal) 
reservoir  shown  in  Figure  5  and  discharged  through  5.076  cm 
(2  in)  diameter  tubing  into  a  cylindrical  flow  settling 


chamber  30.46  cm  (12  in)  in  diameter  and  60.91  cm  (24  in) 


tall  located  within  the  60.91  cm  <60.91  cm  -88.83  cm 


Figure  3.  Probe  and  Traversing  Mechanism 


Jet  flow  pump  reservoir 


Crossflow  and  Jet  Loop  Piping  Diagram 


(24  in  •  24  in  •  35  in)  inlet  chamber  shown  in  Figure  6.  The 
settling  chamber  was  sealed  at  its  bottom  so  that  the  water 
spilled  from  its  top  into  the  inlet  chamber  through  honey¬ 
combed  flow  straighteners  to  reduce  turbulence  and  evenly 
disperse  the  flow.  To  further  reduce  turbulence,  the  flow 
was  broken  by  another  stack  of  honeycombed  flow  straighteners 
and  a  layer  of  fiberglass  filter  material  located  immediately 
above  the  normal  operating  water  level.  The  flow  next 
entered  a  24.4  cm  ■  32.39  cm  ■  182.9  cm  (9.625  in  •'•12.75  in  ■  72  in) 
flume  shown  in  Figure  7  through  a  vertical  section  of  the 
same  honeycombed  material  mentioned  above.  To  avoid  inadver¬ 
tent  spillage  over  the  sides  of  the  flume  during  system 
start-up,  a  5.076  cm  (2  in)  diameter  overflow  pipe  was 
located  in  the  inlet  chamber.  During  normal  operation, 
a  gate  valve  in  this  piping  was  closed.  The  flow  left  the 
flume  through  a  7.614  cm  (3  in)  diameter  pipe  at  its  end 
and  re-entered  the  crossflow  circulation  pump  reservoir.  A 
gate  valve  located  in  this  piping  and  shown  in  Figure  8  was 
used  to  regulate  the  water  level  and  flow  velocity  in  the 
flume.  The  optimum  adjustment  of  this  valve  was  determined 
by  trial  and  error  to  be  closed  two  turns  from  its  fully 
open  position.  Either  a  globe  or  ball  valve  would  have  been 
more  appropriate  for  this  purpose;  however,  neither  was 
readily  available,  so  the  gate  valve  was  used.  The  bracket 
shown  at  the  base  of  the  flume  in  Figure  8  maintained  align¬ 
ment  between  the  flume  and  the  milling  machine.  The  water 
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in  the  flume  was  gradually  heated  by  repetitive  circulation 
through  the  crossflow  pump  and  by  the  addition  of  the  heated 
water  from  the  jet.  To  maintain  a  constant  temperature 
crossflow,  cooling  water  from  a  refrigerated  bath  shown  in 
Fiqure  9  was  circulated  through  a  coil  of  1.269  cm  (.5  in) 
diameter  copper  tubing  located  in  the  flow  settling  chamber. 
Also,  fresh  water  was  added  at  the  flow  settling  chamber 
as  an  equal  amount  was  drained  from  the  crossflow  pump 
reservoir  through  a  1.269  cm  (.5  in)  diameter  pipe.  Cross- 
flow  temperature  was  monitored  by  a  Tvpe-T  thermocouple 
located  in  the  inlet  chamber. 

C.  JET  SYSTEM 

In  reference  to  Figures  4  and  5,  the  jet  flow  pump 
circulated  water  from  a  rectangular  26.27  1  (6.94  gal) 
reservoir  and  discharged  through  1.26  cm  (.5  in)  diameter 
tubing  to  a  33.0  cm  -'50.8  cm  '54.6  cm  (13  in  '<20  in  -  21.5  in) 
head  tank  (Figure  10) .  The  amount  of  flow  to  the  head  tank 
was  regulated  by  a  globe  valve.  Due  to  a  low  flow  rate 
to  the  head  tank,  water  was  also  recirculated  back  to  the 
reservoir  in  order  to  maintain  sufficient  flow  through  the 
jet  pump  to  prevent  overheatinq  it.  A  constant  water  level 
was  maintained  in  the  head  tank  by  a  stand  pipe  which  allowed 
overflow  back  to  the  reservoir.  Sufficient  flow  into  the 
tank  was  maintained  to  make  sure  that  it  slightly  overflowed 
continuously.  Water  drained  from  the  bottom  of  the  head 
tank  through  1.26  cm  (.5  in)  diameter  tubing  and  passed 
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Figure  9.  Refrigerated  Bath 
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through  a  rotometer,  a  .95  cm  (.375  in)  tubing  reducer,  a 
dye  injection  system  (Figure  11) ,  a  water  heater  (Figure  12) 
consisting  of  approximately  6.09  m  (20  ft)  of  .95  cm  (.375 
in)  diameter  copper  tubing  coiled  in  a  heated  bath  and 
finally  a  7.144  cm  (.28125  in)  nozzle  which  discharged  into 
the  bottom  of  the  flume.  Flow  was  controlled  by  pinching 
the  tubing  between  the  heater  and  the  nozzle  with  surgical 
clamps.  Drainage  from  the  crossflow  reservoir  discussed  in 
Section  I1I.B  was  used  to  replenish  the  jet  reservoir.  The 
dye  injection  system,  used  in  photographing  the  jet,  was 
located  approximately  8.23  m  (27  ft)  upstream  of  the  nozzle 
to  minimize  any  disturbance  to  the  jet  that  it  might  have 
caused.  The  majority  of  this  distance  was  taken  up  by  the 
heating  coil  mentioned  above.  The  vertical  distance  between 
the  top  of  the  stand  pipe  in  the  head  tank  and  the  tip  of 
the  nozzle  in  the  flume  was  2.2  m  (86.5  in)  which  equated  to 
21.56  KPa  (3.127  psig) .  Jet  flow  temperature  was  monitored 
by  a  Type-T  thermocouple  located  within  the  jet  flow  tubing 
approximately  1.167  m  (46  in)  from  the  nozzle. 

D.  TEMPERATURE  PROBE 

Measuring  temperatures  in  a  buoyant  jet  with  a  thermo¬ 
couple  is  intrusive.  To  reduce  the  probability  of  distorting 
results,  steps  were  taken  to  minimize  the  cross-sectional 
area  of  the  temperature  measuring  device  as  seen  by  the  flow 
of  the  jet.  A  .02  54  mm  (.001  in)  diameter  Ty:.e-E  micro¬ 
thermocouple  was  selected.  The  suspension  device  for  the 


microtherinocouple  had  to  be  rigid  and  have  a  small  cross- 
sectional  area,  for  reasons  discussed  above,  as  well  as  be 
an  electrical  insulator  to  prevent  interference  with  the 
thermocouple  performance.  A  glass  annulus  approximately 
1.45  mm  (.057  in)  in  diameter  and  11.27  cm  (4.4375  in)  in 
length  was  chosen.  One  lead  of  the  thermocouple  was 
threaded  through  the  annulus  and  the  other  was  glued  with  a 
fast  drying  modelers'  glue  along  the  outer  surface,  allowing 
the  microthermocouple  junction  to  protrude  slightly  from  the 
tip  of  the  annulus.  The  leads  at  the  opposite  end  of  the 
annulus  were  welded  to  .0762  mm  (.003  in)  diameter  wire 
which  subsequently  was  connected  to  28  AWG  extension  wire 
to  the  computer.  The  annulus  was  mounted  as  shown  in  Figure 
13.  Henceforth,  this  device  will  be  referred  to  as  the 
probe . 

E.  PROBE  ACTUATOR  ASSEMBLY 

The  cross  sectional  area  of  the  probe  as  seen  by  the  jet 
was  further  reduced  by  orienting  the  probe  tangentially  to 
the  trajectory  of  the  jet  as  shown  in  Figure  14.  This 
photograph  indicated  that  the  probe  created  no  noticeable 
interference  with  the  jet  hydrodynamics.  Probe  orientation 
was  accomplished  by  the  linkage  assembly  shown  in  Figure  15. 
The  fixed  end  of  the  probe  was  hinged  to  a  streamlined  tube 
23.495  cm  (9.25  in)  long  with  a  maximum  width  and  depth,  as 
seen  by  the  jet,  of  3.175  mm  (.125  in)  and  6.35  mm  (.25  in) 
respectively.  It  was  rigidly  connected  to  the  mounting 
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A.  Drive  Motor  Pulley 

B.  Cable  and  Tumbuckle 

C.  Stroke  Pulley 

D .  Linkage  Rod 

E.  Mounting  Base 

F.  Probe  Arm 

G.  Mounting  Bracket 

H.  1/2-13  Nut 

I.  1/2-13  Hollowed  Stud 

J.  Streamlined  Tube 

K.  Hinge  Assembly 

L.  Probe  (Glass  Annulus) 


bracket  as  shown  in  Figure  15  which  was  connected  to  the 
probe  arm  shown  in  Figure  3  by  a  single  stud  which  allowed 
pivoting  of  the  probe  from  side-to-side .  The  stud  was  also 
hollowed  so  that  a  linkage  rod  could  extend  from  the  hinge 
assembly  through  the  tube  and  stud  to  a  stroking  pulley 
which  was  rotated  by  a  small  motor.  The  hinge  assembly  and 
the  stroking  pulley  were  spring  loaded  to  reduce  hysteresis. 
As  shown  in  Figures  16  and  17,  the  1.5  VDC  motor,  geared  to 
one  rpm,  was  directly  coupled  to  a  potentiometer  as  well  as 
the  drive  pulley.  The  potentiometer  was  configured  in  a 
voltage  divider  such  that  the  amount  of  motor  rotation,  and 
ultimately  the  degree  of  probe  deflection,  was  proportional 
to  the  potential  difference  sensed  across  the  potentiometer. 
Limit  switches  were  installed  at  the  stroke  pulley  as  shown 
in  Figures  16  and  18  to  prevent  damage  to  the  linkage 
assembly  due  to  over-rotation. 

F.  MICROCOMPUTER  INTERFACE 

The  data  collection  process  consisted  of  adjusting  the 
probe  angle  of  deflection,  traversing  the  three-dimensional 
positioning  platform  and  measuring  temperature  profiles. 

All  of  the  mechanisms  which  controlled  these  events  were 
interfaced  to  an  HP-9826  computer  shown  in  Figure  19  through 
an  HP-6942A  multiprogrammer  which  performed  high  speed  analog 
t.o-digital  conversions  and  ultimately  provided  control 
signals  to  govern  relays  within  the  system.  Refer  to  MAIN_T 
in  Appendix  B  for  the  microcomputer  software  which  directed 
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As  discussed  in  Section  III.E,  a  potentiometer  con¬ 
figured  as  a  voltage  divider  provided  probe  angle  feedback 
to  the  computer  as  shown  in  Figure  20.  The  direction  of 
motor  rotation  was  controlled  by  the  computer  through  Type- 
MRR1CDL  replays  connected  as  shown  in  Figure  21.  When  the 
probe  was  at  a  desired  position,  a  5.0  VDC  signal  was  applied 
to  pins  6L  and  6R  which  allowed  both  relays  to  assume  the 
normally  closed  (NC)  positions  which  opened  the  power  circuit 
to  the  motor.  When  it  was  desired  to  rotate  the  motor  clock¬ 
wise,  pin  6R  was  grounded  which  resulted  in  5.0  VDC  applied 
across  the  coil  in  the  right-hand  relay.  This  caused  the 
relay  to  assume  its  normally  open  (NO)  position  resulting  in 
a  1.5  VDC  signal  at  terminal  B  of  the  motor  causing  it  to 
rotate  in  the  clockwise  direction.  The  left-hand  relay  was 
activated  in  a  similar  manner  for  counterclockwise  rotation. 
Refer  to  PROBE_SUBS  in  Appendix  B  for  probe  positioning 
software . 

2 .  3-D  Positioning  Platform  Movement 

Positioning  platform  movement  was  controlled  in  a 
manner  similar  to  the  probe  and  was  discussed  in  detail  by 
Nickodem  [Ref.  2J  .  Refer  to  MTR_SUBS  in  Appendix  B  for 
associated  HP-9826  software. 

3 .  Temperature  Data  Collection 

Three  thermocouples  were  monitored  in  the  data  col¬ 
lection  process.  A  Type-T  thermocouple  located  in  the  inlet 
chamber  measured  the  ambient  fluid  temperature  in  the  flume. 
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a  Type-T  thermcouple  located  in  the  tubing  between  the  heater 
and  the  nozzle  measured  nozzle  temperature  and  a  Type-E 
microthermocouple  in  the  probe  measured  the  temperature  in 
the  jet.  The  EMF's  generated  by  these  thermocouples  were 
amplified  by  "Omega  Omni-Amp  IIB"  millivolt  amplifiers  shown 
in  Figure  22  prior  to  entering  the  multiprogrammer  for 
ana  log- to-digi tal  conversion  and  eventual  transformation  to 
temperature  readings.  Fourth-order  least  squares  coeffi¬ 
cients  for  this  conversion  were  taken  from  Beckwith  [Ref.  5] . 
Operation  of  the  crossflow  circulation  pump  created  suffi¬ 
cient  electrical  interference  to  distort  the  thermocouple 
signals.  This  problem  was  corrected  by  applying  a  thin 
coating  of  silicon  sealant  to  the  Type-T  thermocouple  junc¬ 
tions  and  by  connecting  the  crossflow  circulation  pump 
casing,  the  nozzle  and  the  jet  tubing  in  the  vicinity  of  the 
nozzle  thermocouple  to  a  common  ground.  Because  the  jet 
tubing  was  plastic,  it  was  necessary  to  manufacture  a  brass 
"T"  connector  as  shown  in  Figure  23  which  was  grounded  and 
located  in  close  proximity  to  the  thermocouple  junction. 

Refer  to  T_SUBS  in  Appendix  B  for  associated  software. 
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Figure  22. 


Thermocouple  Amplifiers  and  HP-6942A 
Multiprogrammer 
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Figure  23.  Nozzle  Thermocouple  Grounding  Arrangement 


42 


IV. 


EXPERIMENTAL  PROCEDURES 


A.  CALIBRATION 

Before  the  data  collection  process  could  begin,  the 
rotometer,  the  thermocouples,  the  probe  and  the  positioning 
platform  had  to  be  calibrated.  The  detailed  steps  taken 
are  discussed  below. 

1 .  Rotometer 

With  a  constant  level  maintained  in  the  head  tank 
and  the  jet  tubing  disconnected  from  the  nozzle  and  elevated 
to  the  same  height  as  the  top  of  the  nozzle,  five  100  ml 
samples  were  drawn  through  the  rotometer  and  timed  to  the 
nearest  0.01  second  at  each  rotometer  reading  from  10° 
to  751  in  5%  increments.  Flowrates  and  standard  deviations 
in  ml/s  are  shown  in  Table  1. 

2 .  Thermocouples 

Since  nozzle  and  probe  temperatures  were  to  be 
normalized  by  the  ambient  temperature,  the  nozzle  and  probe 
thermocouples  were  calibrated  relative  to  the  ambient  thermo¬ 
couple  by  using  the  microcomputer  program  T_CAL  in  Appendix 

B.  The  procedure  followed  is  outlined  in  the  initial  com¬ 
ments  of  the  program.  Coefficients  for  first  order  curve 
fits  were  solved  by  the  least  squares  method  with  the 
mainframe  programs  TCAL  and  TFIT  found  in  Appendix  C. 
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Probe 


As  the  probe  assembly  was  being  developed,  it  was 
convenient  to  test  its  suitability  with  the  probe  calibration 
panel  shown  in  Figure  24.  Resistance  changes  across  the 
potentiometer  were  recorded  for  varying  degrees  of  deflec¬ 
tion.  Analysis  of  this  information  led  to  improved  designs 
from  the  standpoint  of  reduced  hysteresis  and  repeatability. 
The  microcomputer  program  PROBE_CAL  in  Appendix  B  was 
developed  to  enable  calibration  of  the  final  design  after 
it  was  installed  in  the  system  as  shown  in  Figures  25  and 
26.  The  calibration  procedure  is  outlined  in  the  preliminary 
comments  of  the  program. 

4 .  3-D  Positioning  Platform 

The  positioning  platform  was  calibrated  in  a  manner 
that  placed  the  tip  of  the  probe  at  desired  locations  within 
the  flume  relative  to  the  tip  of  the  nozzle.  Referring  to 
the  coordinate  system  illustrated  in  Figure  2,  the  center 
of  the  nozzle  was  defined  as  (0,0,0)  in  xy z-coordinates . 

The  following  relationships  apply  to  the  probe  geometry 
shown  in  Figure  27: 

X(real)  =  X  -  R  cos  t  tan(*r/4  -  a/2) 

O  <3 

Y(real)  =  Y  +  r  (1  -  cos  .)  +  R  cos  a 
op  a 


Z  ( real )  =  Z  -  r  sin  , 

o  p 
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Figure  27.  Probe  Positioning  Geometry 
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(X  ,Y  ,2  )  was  the  position  of  the  millina  machine  bed 
o  o  o 

when  the  tip  of  the  probe  was  at  position  (0,0,0)  with 

:  =  0  and  t  =  .  =  90  decrees.  For  calibration,  the  probe 

and  probe  arm  were  configured  with  these  settings  as  shown 

in  Figure  28.  With  measured  values  of  r  and  R  entered 

pa 

into  the  microcomputer  program  MAIN_T,  the  calibration  was 
accomplished  by  the  program  MOTOR_CAL.  The  step-by-step 
procedure  followed  was  outlined  in  the  subprogram  "SUB 
Calibrate."  As  illustrated  in  I igures  3,  27  and  28,  the 
length  of  could  be  modified  to  compensate  for  adjustments 

of  a  and  3  to  positions  other  than  90  degrees.  Decreasing 
x  increased  the  distance  along  the  Y-axis  in  which  the 
probe  could  be  positioned.  Increments  of  x  and  were 
scribed  on  the  top  of  the  base  and  at  the  tip  of  the  probe 
a^m  in  Figure  28  to  accommodate  this  change,  if  desired. 

The  program  MAIN_T  queried  the  user  for  the  value  of  \  and 
assumed  t  =  2.  The  calibration  software  also  established 
position  limits  to  prevent  driving  the  probe  into  the  sides 
of  the  flume. 

B.  PRELIMINARIES 

Crossflow  velocity  was  determined  by  injecting  blue  food 
coloring  into  the  flow  and  timing  its  travel  through  a  1.0m 
interval.  The  average  of  several  trials  indicated  the 
velocity  was  .1)0  m/s  (.427  ft/sec)  with  the  flume  outlet 


valve  closed  two  turns  from  its  fully  open  position. 


Photographs  of  jet  profiles  as  shown  in  Figure  29  we  re¬ 
taken  to  determine  jet  trajectory  and  halfwidths  alone  cr. 
streamwise  axis.  This  was  done  by  injecting  blue  feed 
coloring  into  the  jet  flow  as  discussed  in  Section 
The  specific  gravity  of  the  food  coloring  was  found  to  :.o 
considerably  less  than  that  of  water.  To  eliminate 
added  buoyant  effect  this  would  have  had  on  the  jet,  a  smell 
quantity  of  alcohol  was  mixed  with  the  food  coloring  as 
suggested  by  Merzkirch  [Ref.  6] .  The  amount  of  alcohol  added 
was  determined  by  trial  and  error.  As  small  quantities  were 
added  and  mixed,  samples  were  gently  placed  on  the  surface 
of  a  beaker  of  water.  Pure  food  coloring  laid  on  the  sur¬ 
face  and  very  slowly  mixed  with  the  water.  As  alcohol  was 
added,  this  buoyant  effect  grew  progressively  less  and  the 
mixture  would  settle  into  the  water.  The  mixture  was  con¬ 
sidered  satisfactory  when  it  no  longer  laid  on  the  surface, 
but  settled  to  some  equilibrium  position  in  the  beaker. 

Slide  photographs  of  the  jets  were  projected  onto  large 
sheets  of  3.175  mm  (.125  in)  grid  graph  paper  and  digitized 
along  approximate  streamwise  axes  and  half-width  trajectories 
A  scaling  factor  was  determined  by  equating  the  projected 
width  of  the  nozzle  to  its  known  outer  diameter  of  7.9375  mm 
(.3125  in).  The  above  data  was  fit  to  the  Michael is-Menter 
Equation  [u-:f.  7]  shown  below  by  the  least  squares  method 
with  the  mainframe  program  JETCURV  in  Appendix  C: 
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Correlation  coefficients  close  to  1  were  consistently 
obtained.  To  determine  positions  within  the  jet  at  which 
to  make  temperature  measurements,  five  evenly  spaced  posi¬ 
tions  per  jet  flow  rate  were  selected  along  the  streamwise 
axis  in  the  zone  of  established  flow.  Data  planes  slightly 
larger  than  the  jet  width  were  centered  at  these  points  and 
oriented  perpendicular  to  the  streamwise  axis.  One  hundred 
data  points  were  selected  in  a  symmetric  square  matrix  with 
points  most  densely  populated  near  the  center.  The  planes 
were  identified  alphabetically  and  in  consecutive  order  from 
"A"  to  "F",  where  "A"  represented  the  plane  nearest  the 
nozzle.  The  positions  were  entered  into  the  microcomputer 
and  stored  by  plane  on  a  floppy  disk  by  the  program  LOAD_XYZ 
in  Appendix  B.  Accompanying  each  data  point  was  a  probe 
deflection  angle  used  to  orient  the  probe  parallel  to  the 
path  of  the  jet  to  minimize  interference.  This  angle  was 
determined  by  evaluating  the  first  derivatives  of  the 
equations  developed  for  the  streamwise  axis  and  the  half¬ 
widths  and  performing  an  interpolation  based  upon  the  data 
points'  position  relative  to  the  two  curves. 

C.  DATA  ACQUISITION 

The  flow  systems  were  placed  into  operation  and  the 
ambient  and  nozzle  temperatures  were  monitored  with  the 
program  T_SCBS  to  evaluate  system  stability  and  readiness 
for  data  acquisition.  The  system  usually  took  approximately 
two  hours  to  come 
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into  ecu i 1 ibr ium . 


This  could  be  monitored 
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by  watching  the  jet  nozzle  temperature.  When  conditions 
were  stable,  data  consisting  of  two  hundred  probe,  ten 
ambient  and  ten  nozzle  temperature  samples  per  position  was 
collected,  one  plane  at  a  time,  by  the  program  MAIN  T. 

The  following  information  was  stored  on  a  floppy  disk  for 
each  data  point:  x,  y,  and  z  coordinates;  mean  probe, 
nozzle  and  ambient  temperatures  and  the  standard  deviation 
of  the  probe  measurements.  The  data  was  transferred  to 
the  mainframe  computer  by  using  a  modem,  the  microcomputer 
program  SEND_DATA  and  the  mainframe  program  GRAB. 

D.  DATA  REDUCTION 

The  raw  data  was  organized  into  a  more  usable  format 
by  the  mainframe  program  TDATA  which  also  converted  the  XYZ 
coordinates  into  the  XSW  system  shown  in  Figure  2.  The 
resulting  data  was  selectively  sent  to  the  program  C0NT0UR4 
which  applied  calibration  coefficients  to  the  temperature 
data  and  normalized  it  in  the  following  manner: 


n  a 


where  T  was  the  jet  temperature  as  measured  bv  the  probe, 

p  J  1 

T  was  the  ambient  fluid  temperature  and  T  was  the 
a  ^  n 

temperature  of  the  jet  within  the  nozzle. 

Contour  plots  of  this  ir formation,  generated  by  the 
CONTOUR  option  of  the  graphics  package  DISSPLA  [Ref.  8], 
are  presented  in  Figures  30  through  35. 
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Figure  32 


Plane  B  Temoerature  Contour  Plot 
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Figure  34.  Plane  D  Temperature  Contour  Plot 
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RESULTS 


E  . 

The  contour  plots  found  in  Figures  30  through  35  support 
the  conjectures  in  Section  II. C  concerning  the  effects  of 
the  crossflow.  Figure  30,  a  plot  of  plane  A  centered  on 
the  streamwise  axis  46  degrees  from  horizontal  and  7.327  mm 
(.288  in)  downstream  of  the  nozzle  shows  significant  distor¬ 
tion.  Figures  31  through  35  show  planes  A  through  E  sequen¬ 
tially  plotted  on  the  same  scale  in  order  to  observe  overall 
jet  behavior.  Plane  E  was  located  87.313  mm  (3.438  in) 
downstream  of  the  nozzle  and  86  degrees  from  horizontal. 

It  can  be  seen  that  as  the  jet  traveled  further  downstream, 
the  distorting  effect  grew  progressively  less,  as  expected. 

The  rate  of  heat  transfer  from  the  jet  to  the  ambient 
was  calculated  for  each  plane  utilizing  the  temperature 
distribution  matrix  generated  in  the  program  CONTOUR4  in 
Appendix  C  and  the  following  relationship: 


Q 


m  n 

l  l 


i=l  j=l 


U 


c  A. .  U. 

P  13  13 


T.  . 
13 


where  o  was  the  relative  density  of  the  jet,  A. .  was  the  area 

13 


of  each  matrix  segment,  IF  ^  was  the  velocity  in  each  segment, 


c  was  the  specific  heat  of  the  jet  and  T. .  was  the  tempera- 
ture  in  each  segment.  Velocity  was  measured  along  the 
streamwise  axis  by  a  laser  Doppler  velocimeter  (LDV).  It 
appeared  to  be  nearly  constant  in  the  region  of  the  jet 
observed.  The  mean  and  standard  deviation  was  44.875  mm/s 


and  1.8  mm,  s  respectively.  By  assuming  the  Gaussian  profile 
shown  in  Section  II. B,  velocity  was  determined  at  radii 
corresponding  to  each  segment  in  the  matrix  mentioned 
above.  The  rates  of  heat  transfer  are  shown  in  Table  2. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  objective  of  this  thesis  was  to  develop  a  computer- 
aided  data  acquisition  system  and  construct  a  microthermo¬ 
couple  probe  to  be  used  by  follow-on  students  to  study 
temperature  distributions  in  turbulent  buoyant  jets.  Sample 
data  was  taken  to  verify  system  operability.  Based  on 
results,  the  system  performed  in  a  satisfactory  manner  and 
will  be  an  invaluable  tool  for  subsequent  studies  of  buoyant 
jets  in  a  crossflowing  ambient. 

Data  point  positions  were  hand  calculated  and  loaded 
into  the  microcomputer  with  the  program  LOAD_XYZ.  This  was 
an  extremely  time  consuming  task  and  distracted  the  user 
from  defining  more  than  100  data  points  per  plane.  The 
system  can  be  greatly  improved  with  the  addition  of  micro¬ 
computer  software  that  would  automatically  determine  and 
load  data  point  positions.  Data  point  population  could  then 
be  increased  with  ease  which  should  result  in  smoother  con¬ 
tour  plots  and  more  accurate  heat  transfer  calculations. 

It  was  necessary  to  continuously  add  fresh  water  to  the 
crossflow  system  as  an  equal  amount  was  drained  from  it  in 
order  to  maintain  the  crossflowing  ambient  at  a  constant 
temperature.  This  was  because  the  cooling  coil  located 
within  the  flow  settling  chamber  was  inadequate  to  compen¬ 
sate  for  the  heat  added  to  the  system  by  the  jet  and  the 


crossflow  circulation  pump  by  itself.  Although  it  was 
possible  to  maintain  the  ambient  constant  within  1.4  C  b 
this  method,  in  the  interest  of  conserving  water,  it  is 
recommended  that  either  a  larger  capacity  chilled  water 
bath  or  a  cooling  system  that  circulates  a  refrigerant 
rather  than  water  be  appropriated  for  this  purpose. 
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APPENDIX  A 


UNCERTAINTY  ANALYSIS 

Experimental  uncertainty  was  analyzed  in  accordance  with 
the  guidelines  set  forth  by  Holman  [Ref.  9] .  Uncertainties 
in  the  primary  measurements,  based  upon  manufacturer  speci¬ 
fications  and/or  the  number  of  significant  digits  which 
could  be  read,  follow: 

1.  Time  (rotometer  calibration):  .01  s 

2.  Volume  (rotometer  calibration):  1.0  ml 

3.  Rotometer  reading:  1.0  percent 

4.  3-D  positioning  platform  resolution:  .1524  mm 

5.  Probe  deflection  angle  resolution:  .5  degrees 

6.  Thermocouple  resolution:  .0099  °C 

7.  Thermocouple  time  constant:  .004  s 

8.  Time  (crossflow  velocity):  .01  s 

9.  Length  (crossflow  velocity)  :  .05  m 

Based  on  the  above,  the  uncertainty  of  the  flow  rate  of  the 
jet  was  estimated  to  be  1.0  ml/s.  Maximum  uncertainty  in 
the  position  of  the  probe's  tip  due  to  the  uncertainty  in 
deflection  angle  was  determined  to  be  .983  mm  (.0387  in). 
Combined  with  the  resolution  of  the  positioning  platform, 
the  tip  of  the  probe  was  positioned  with  an  uncertainty  of 
1.135  mm  (.0447  in)  in  each  plane,  or  with  an  overall 
uncertainty  of  1.605  mm  (.063  in)  in  three-dimensional 
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APPENDIX  B 

MICROCOMPUTER  PROGRAMS 

mhin  T 


■i 


;o 

•  'Si; 

3  r 09 ran  coordinates  tne 

ui; 

'  data-* 

aking  evolution  tor  op-isur  . 

- 

50 

'  temperature  JistriDutions  in  ouovant 

m. 

60 
-  *1 

*  jets. 

i 

t  • j 

00 

i 

30 

!  *  1  . 

Load  all  SuDproqruns 

" 

100 

»  2 . 

Input  oesired  positions  t run  i 

!  '0 

t 

disk  file  of  the  form  ""UNtX” 

120 

i 

lvalue  "XX"  is  tne  run  niinDer> 

1  20 

f  3. 

■Move  'he  3_D  positiomnq 

• 

i  40 

\ 

platform  to  each  Position 

150 

!  4 . 

Aiiqn  tne  thermocouple  prone  Jith 

1 60 

i 

the  jet  s t reami)  1  se  axil 

1/0 

I  5. 

Obtain  200  temperatures  at  each 

130 

i 

position  ana  compute  tne  mean  ana 

130 

i 

standard  deviation  isd) 

200 

»  6. 

Write  to  disn: 

2 ' 0 

» 

a .  T  ( C ) 

• 

220 

t 

b .  sd 

230 

» 

C.  x  .  /  .  t  in  (mm) 

240 

» 

not  tie  centerline  13  <0.v.v> 

250 

» 

250 

!  Probe 

and  arm  dimensions  in  inches: 

2/0 

?  a . 

Shortest  arm  lenuth  ■  .A  . .. 

230 

!  b. 

Lonqest  arm  lenqth  -  24.3 

230 

i 

• 

200 

Lenqth_proDe*4 .4375 

310 

320 

Leng  th_arm 
♦ 

-20.0 

330 

OPTION  BASE!  1 

340 

DIM  CoeT  (  >2) .X<50G> , T  <  500  > ,  Z <500 > .  P  rot>e_anq  i  e  <500  > 

350 

LOADSUB  ALL  PROM  "T  GOBS" 

1 

360 

LOADSUB  ALL  FROM  "MTR  SUBS" 

- 

370 

LOADSUB  ALL  from  "PROBEJSUCS" 

•1 

280 

CALL  Retrieve  coef <Coef <* > , "motor_coet " ) 

330 

t 

4[J0 

!  2.  Input  desired  positions  from  disk 

410 

1 

■ 

420 

1=1 

- 

430 

INPUT 

"Angle  of  ar'n  relative  To  ♦Y-iXis?" 
"F  1  ienaine  for  dos  i  t  4 on  i  ng  daUv'  ,r  i 

.Angle  arn 

1 

4u0 

INPUT 

1  pTfcime  i 

■ 

450 

ASSIGN 

*F  i  le4  TQ  Filenames 

460 

Go_on  : 

♦ 

4  70 

ENTER  ?F 1 1 e4 ; X < I ) , Y  < ; >  .  Z <  I  > 

480 

1 

4  30 

IF  ’<(1)0-100  THFN 

500 

IF  x<  D--338  r MEM 

510 

P  rode  _anq i e  < I ) «  Y <  I ) 

520 

D  anq ie*T  <  i  > 

530 

ELSE 

• 

540 

Probe  anqie(I)=P  anqie 

550 

I  •  I  ♦  1 

OmD 
5  70 

J  r‘. ') 

END  Ir 

* 

GOTH  ( 3 0  on 

5  20 

end  if 

.0  0 
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650 

608 

670 

-S30 

630 

700 

710 

720 

730 

740 

750 

760 

770 

760 

790 

300 

310 

820 

330 

e_anqle< 
840 
850 
860 
870 
380 
8 '30 
900 
910 
920 
930 
940 
350 
960 
370 
980 
990 
1000 
1010 
1  020 
1030 
1040 
1050 
1060 
1070 

1  n«o 
1090 


bu  .  0  bo  _o n 
t '< D  IF 

t 

ASSIGN  «File4  TO  * 

N  i  terns  -  I  -  1 

i 

'3.  Begin  loop  to  tax  e  data  at  each  point 

i 

INPUT  "NAME  OF  FILE  WH ERE  DATA  IS  TO  BE  STORED’" , F i 1 enaroe 1 $ 
Records*  <  N  i  terns *8*  7/256 )  *2 
CREATE  BOAT  F i 1 ename 1 $, Records 
ASSIGN  4F  i  lei  TO  FilenamelS 


FOR  I_position“1  TO  Nitems 

i 

'a.  Move  milling  machine  and  move  the 
'  thermocouple  probe. 

i 

I ■ I _pos l t ion 

CALL  Move_ldv_to<X(I) ,T(I) ,Z<I) ,Lenght_arm, Leng t h_p robe , Ang 1 e_arm 
I )  ,Coef  < *  > ) 


!b.  Obtain  mean  temperature  and  sd 

i 

call  T_coup le<  T , "PROBE" , "C" . 200 . S  t_dev  > 

i 

!c.  Measure  ambient  and  jet  temps 

t 

CALL  T._couple<  T_amb lent," AMBIENT ","C", 10.Sd> 
CALL  T_coupie(T_no2zle," NOZZLE ","C".10,Sd> 

i 

!d.  Write  ail  information  to  disk 

i 

OUTPUT  SFilel ; X <  I ) , Y (  I ) , Z < I ) 

OUTPUT  '«F  1 1  e  1  :  T  ,  T_amb  ient,T_nozzle 

output  ?f i i e i : s t_dev 

NEXT  I_pos i t ion 
♦ 

!  4  .  Close  files 

i 

OUTPUT  <>Fi  lei  ;  -  1  00 
ASSIGN  *Fi  lei  TO  • 


'5.  Terminate  program 

t 

PRINT  "All  done!" 


PROBE  SUBS 


: ) 

40 
50 
50 
70 
50 
00 
'  00 
!  10 
120 
no 

140 

no 
160 
:  70 

no 
no 
200 
210 
220 
230 
240 
250 
260 
270 
280 
230 
300 
310 
320 
330 
340 
350 
360 
2  70 
380 
330 
400 

410 

420 
4  30 
440 
4  50 
4b  0 
470 
480 

4  90 
500 
510 
520 
5.30 
Sul 
550 
550 

5  70 


'  '.H5  program  moves  trie  temperature  oroDe 

'  to  desired  angies  of  oe  r  i  ec  t .  or, . 

t 

'  NOTE:  Calinration  roe f ( i c  i  en t s  are 
!  entered  i n  303  Read_anqie. 

!  beginnmq  at  line  290. 

i 

COLL  Read _ang 1 e< Ang le > 

oRINT  "The  prone  is  presently  at "  :  Ang  i  e  ;  "  degrees  from  nonzonta 
INPUT  "rinar.  is  your  desired  angie  *or  the  p  r  oDe  ,  5es  i  r  ec  _ang  i  e 
CAl  L  3 rooe _move < Des i red _anq i e > 

go t o  no 

EMI) 

t 

i 


SUB  Read_anq 1 e< Ac t ual _ang 1 e > 

,This  program  reads  tne  present  angle 
lot  tne  proDe  and  returns  it 
'  0  degrees  *  Pot iiontai 
190  degrees  *  vertical,  dout. Liard 


1A0,  At  and  A2  -  coefficients  for  a 
Isecond  order  curve  fit  of  ml/  vs  angie 
!of  deflection  data. 

i 

1  Coe f f ic len ts  for  6  dune  data  follow. 


A0--23. 3657824056449 
A1 *3 . 099955062331 b4E- 2 
A2= -  1 .229223054561 J9E -6 
'  !  F0RHA T  A/D  CARD 

'OUTPUT  723  :  "CC  ,  I  f”  IClear  Relay  Card 
OUTPUT  723:  "CC  .  3 T"  'Clear  A/D  '.arc 
OUTPUT  723 :  "CC .  7T"  'Clear  Digital  Card 
OUTPUT  723;"5F ,3.8.3. 1 .25.12T" 

OUTPUT  723 : "DB , 1 , 1 0 . 1 T"  '  CLOSE  RELAT 
OUTPUT  723  :  " TP . 3T"  1  START  9/D 

ENTER  7230  I  :  0 

Actual  ang 1 e* A0 ♦ A 1 “V *A2*V*V 
SUBEND 


SUB  Prone_move( Des l red _anq 1 e > 

T  ft  i  s  suDProgram  moves  the  orotie 
desired  angle 


to  the 


0.  Check  to  see  if  the  anqie  is  m 
an  acceptao.e  range. 


OUTPUT  723 : "OP .1,0 T"  'Clear  Relays 
CAll  C 1  ear _sc f een 


IF  Des  i  red  ang  l  e>4|)  THEM 


r>  0  0 

366?  3400.1 

-:;o 

-‘’In'  "Desired  angie  e<ceeds  30  degrees 

“J  w  o 

SU6E  <  I  ’ 

n  :>  0 

END  IF 

1 

^  ’j 

350 

IF  Desired  angle<0  THEN 

550 

BEE?  3400.’ 

570 

3EEP  3800.1 

o80 

BEEP  3600.1 

630 

P P I N T  "Desired  angle  is  negat l ve • 1 ! " 

7  00 

SUBEXIT 

710 

END  IF 

720  ’ 

7  30  ’ 

l.a,  Clear  the  digital  output  card. 

7  4  0  ! 

b.  Format  the  A/D  card. 

750  ’ 

c.  Close  the  relay  that  corrects 

760  ! 

the  probe  potentiometer  to  the 

7  70  • 

A/D  conve rter . 

780  ’ 

790 

OUTPUT  723;”OP,7,0T” 

300 

OUTPUT  723 : "SF , 3 , 3 . 3 , 1 .25, 3T" 

310 

OUTPUT  723:"CC, 1 T" 

320 

OUTPUT  723;  "OB, 1 .10,11" 

330  ’ 

840  ’ 

2.  Define  the  acceptable  tolerance 

850  ! 

in  the  angle  (degrees). 

360  ’ 

3  70 

Tol erance- . 5 

380  ’ 

390  ’ 

3.  Control  loop. 

900  ’ 

910  Repeat: 

CALL  Read  angleCActual  angle) 

920 

PRINT  "ANGLE  *" ; 

930 

PRINT  USING  "DDD . DD” ; Ac tual _ang 1 e 

940 

BEEP  Ac t ual_ang i e*  1  00 ,  .  05 

950 

Angle  e r ro r - < Oes i red  angle-Actual  angle) 

1 

960 

3  70 

IF  A8S < Ang le_e r ror >> To  1 erance  THEN 

i 

980 

390 

IF  Ang 1 e_e r ro r > *0  THEN 

1000 

Directions-" Down" 

1010 

ELSE 

1020 

Directions- "Up " 

1030 

END  IF 

10  40 

1 

1050 

CALL  No tor_go < D i rec t i onS ) 

1  060 

t 

’070 

GOTO  Repeat 

1080 

END  IF 

1030  OUTP’JT 

723 ; "OP , 7 , 0T" 

’’00  QlJfp"'  723  ;  "CC  ,  ■  T"  ’Clear  relay  card. 
1110  SiJBE MD 

f 

1130  ’ 

11^0  1 

1150  S1JB  No  t  o  r_3o  <  D  t  r »c  t 1  onS ) 


'200  IF  D  1  rec  t  ionij‘'‘Doun"  ~HEN 

!2’0  LDit*7 

1220  END  IF 

1230 

l  240  OUTPUT  723 ; "OP , 7" :2  Lb  it;" 
1250  SUBEND 
1260  ' 

1270  ' 

1280  * 

1230  SUB  Clear_screen 
1300  1  Clear  the  CRT. 

1310  1 

1320  OUTPUT  2  USING  "*  ,  B‘* : 255 , 75 

1330  GCLlAk 

1340  SUBEND 
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MTR  SUBS 


!  me  ^allowing  series  of  suoroutines 
'  ire  utilized  to  caiiDrate  me  Positioning 
'  iidt'orn  ana  ultimately  to  move  me  orooe 
1  tip  to  desired  positions  .jit’iin  tn^  jet. 


SO 

SUB  D raw_f 1 ume 

■00 
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!  ’  n 

f  Draw  the  Buoyant  Jet 
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■  20 
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I  60 
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i  ;o 
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Ml) 

1  Draw  the  too  view 

'  30 

IBP AW  26.0 

200 

I DRAW  0  ,  “  n 

210 

I DR AW  -36,0 

220 

I  DRAW  0  .  - ' 0 

230 

i 

2«0 

move  6.o 

2S0 

I D  R  a  W  0  .  '  0 

260 

MIME  24.0 

270 

IDRAW  0.  -10 

230 

'  iabel  "Top1'. 

2°0 

hOVE  0.10 

200 

:s:ze  5 

310 

LABEL  ''Top” 

3  20 

'  Draw  the  side  view. 

2  30 

MOVE  6.22 

.24  0 

I  DRAW  36 , 0 

3S0 

IDRAW  0.13 

360 

IDRAW  -26.0 

2  70 

IDRaw  0,-13 

330 

MURE  6.0 

330 

IDRAW  0  .  '  3 

a  00 

MURE  24,0 

4  10 

IDRAW  0.-13 

4  J  0 

'  .aDel  "Side" 

4  20 

2SIZE  5 

4  40 

'•ORE  0 . 28 

4  SO 

..ABEL  -S.de" 

4m) 

\ 

470 

1  Label  the  picture. 

46  U 

» 

4  60 

'WE  M.35 

Slid 

f '  r*  r  7  £•  j 

310 

.[ABEL  "BUOYANT  ^ET  E 

=  Mi 

i 

3  30 

’  Eit  on  the  no  1 1 1 e . 

3  41) 

f 

350 

'•TIRE  14.22 

5  r>  6 

IDRAW  0.2 

5  70 

IDRAW  .25.0 

5  >  >  1 1 

ir-W-w  1.2 

3)0 

1 

-,0D 

:*4nvF  o.-j 

'  n 

•/;  r  ;r 
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550 

,  660 

I  670 

1  630 

530 
700 
710 
720 
730 
740 
750 
760 
770 
780 
730 
800 
310 
320 
830 
840 
350 
360 
370 
880 
890 
900 
310 
820 
330 
940 


CSIZE  3 
LABEL  "o" 

!  Indicate  the  direction  of  flow. 

i 

MOOE  8,11 
I  DR AW  2.0 
IDRAU  -.5. -.5 
[DRAW  .5, .5 
IDRAU  -  .5,  .5 

i 

MOVE  8.27 
IDRAU  2.0 
IDRAU  -.5, -.5 
IDRAW  .5  .  .5 
IDRAU  -  ,5.  .5 

i 

SUBEND 

i 


i 

> 


t 

SUB  Cal ibratetFi lenameS) 

i 

DPT  ION  BASE  1 
DIM  Coef (12) 

i 

f  Calibrate  the  positioners  on  the 


350 
360 
370 
880 
890 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
10  70 
1080 
'090 
1100 
1  n  0 

1  130 
1  130 
1 140 
1150 
1  160 
1  1  70 
1 '30 
1  1  90 
1200 
, ,  0 


milling  machine  movement. 

Onto  disk,  write  out  the  calibration 
coefficients  and  the  hard  boundaries 
that  must  be  observed! 

This  file  will  be  called  "Motor^coef' 


L  . 


Position  the  probe  /olume  at  the 
mail  of  the  tip  of  the  nozzie.  This 
position  is  (0,0,0).  All  readings 
will  be  in  inches.  Read  all  three 
potentiometers.  Ask  the  user  for  the 
notile  outer  diameter  and  compute  the 
zero  position.  Ask  the  user  'or  the 
milling  machine  readings. 


Next,  move  the  bed  to  some  new  posi¬ 
tion  us  i  nq  the  override  switches. 
rak e  readings  from  the  pots  and  ask 
for  the  milling  machine  readings. 
Compute  the  calibration  coefficients. 


Nove  the  bed  to  each  of  the  extremes 
in  the  X.  Y.  and  Z  directions  usmq 
the  override  switches  and  have  the 
user  tell  the  computer  wnen  each  of 
these  boundaries  are  hit.  Enter  each 
of  tnese  onto  the  disk  file. 


D.  Disk  fiie  "k|otor.coei 


1220 
1230 
i  240 


:250 

!  1  ■  x_zero , 

l ope 

’,260 

1  2  .  y  _  z  e  r  o  . 

y  1 ooe 

1270 

1  3 .  i.iero , 

z  _5  iope 

1  2S0 

1  4 .  x _n in,  x 

jnax 

1230 

1  5 .  y_m  in,  y 

1300 

1  b .  z  _m i n ,  i 

^ax 

1310 

i 

1320 

BEEP 

1330 

PRINTER  IS  1 

1340 

gclear 

1350 

OUTPUT  2  USING 

“*  ,8‘*.255 ,75 

1360 

i 

1370 

PRINT  "I.  POTENTIOMETER  CALIBRATION: 

" 

1380 

PRINT  " 

" 

1330 

PRINT  •*  NOTE:  1. 

Probe  must  be  horizontal" 

1  400 

PRINT  "  2. 

Arm  must  De  parallel 

to 

1410 

PRINT  ” 

tne  bed  axis 

1420 

PRINT  ”  3. 

ALPHA  »  BETA 

1430 

PRINT  ” 

1440 

PR INT  "  A ,  Us  inq 

the  override  switches. 

1450 

PRINT  "  position  the  probe  volume  at 

1460 

PRINT  "  the  outer  wall  of  the  tip  of 

1470 

PRINT  ”  the  nozzie. 

1480 

PRINT  ” 

1490 

PRINT  ”  B.  I  will 

need  the  nozzle  O.D. 

and** 

1500 

PRINT”  the  «i 

llmq  machine  positon. 

1510' 

PRINT  ” 

1520 

PR  IN .  ”1  Hit  <con t  >  i 

1530 

PAUSE 

1540 

BEEP  1500  ,  .1 

1550 

INPUT  ”1 .  Nozzie  C 

.0.  ( inches) 7” . No z z le_ 

od 

1560 

BEE3  2000  ,  .  1 

15  70 

INPUT  ”2  .  X  (  m) . 

l ♦  into  f  1  une i ?" , X_! 

1560 

BEE?  2500  .  .  ’ 

1590 

INPUT  ”3.  y  <  in)  , 

i*  alonq  flume  to  the 

right)7 

1600 

BEER  3000  .  .  ’ 

1610 

INP'.J T  "4  .  z  tin). 

!  *  upward!  ,  Z_1 

1620 

i 

1630 

CALL  Read _po t < "X" 

7x  1  ) 

1640 

CALL  Read_po  t  < 

J y  > 

’650  CALL  ^ead _po t ( "2" . 1 l _ ’ ) 

1660  ' 

1670  * 

1680  CALL  C  iear  _sc  reen 

1690  PRINT  ”f. .  Move  the  ml  ling  machine  to  a  new" 

'700  PRIM!  "  position  in  3-D,  by  at  least  6 
1710  PRINT  ”  inches  in  each  direction. 

1720  BEEP  3200. . 1 

1730  INPUT  "X.  Y.  Z  in  i nches?" , X_2 . Y_2 , Z_2 
17  40  COLL  Read _po t ( "X" , V x _2 > 

1  750  COLL  Read, pot <"Y",Vy_2> 

1760  CALL  Read,po t < ”Z" , V z  2) 

1  770  ! 

1780  ’  C.  Calculate  the  calibration  coefficients 

1  730  ! 

i  300  X  jero--t  ( No  z  z 1  e_od/2)  <-7*  J  •<  <  X  _2  - 1  >  /  ( Vx  _2  -Vx  J  ) )  > 

1  ‘  1  ft  X  _3  i  ope  *  <  X_2  -  X  _1  >  /  <  V *  _ J  -  V  »  _  '  ) 
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1820  ' 

1330  Y  zero*-Vy_! *< Y_2-Y_1 >/( Vy_2-Vy  1  > 

1  840  Y  slope*(Y_2~Y_1  )  /  <  Vy_2 -Vy  > 

1350  !  ~ 

1860  Z_zero*-Uz_1*<Z_2-Z_1 >/(Uz_2-Vz_l > 

1870  Z_siope-<Z_2-Z_1)/<Vz_2-Vz_l  ) 

1880  * 

1890  1 

1900  !  D.  Find  the  physical  boundaries  for  each 

1910  !  d  i  rec  t ion , 

1920  ! 

1930  CALL  Clear_screen 

1940  PRINT  "1,  Move  the  millinq  machine  to  the 
1950  PRINT  "  minimum  value  of  x  .  Hit  <cont> 
1960  PRINT  "(Away  from  the  flume,  backwards; 

1970  PAUSE 

1980  CALL  Read_po  t  < "X” , V ) 

1990  X_m  i  n=-X_ze  ro'*'X_s  iope*U 

2000  PRINT  "2.  Move  the  milling  machine  to  the 

2010  PRINT  "  maximum  value  of  'x'  .  (Towards 

2020  PRINT  ”  the  flume).  Hit  <cont>  . 

2030  PAUSE 

2040  CALL  Read_po t < "X" ,  V) 

2050  X _max *X _ze ro *X_s 1 ope*U 

2060  PRINT  "3.  Move  the  milling  machine  to  the 

2070  PRINT  "  minimum  value  of  'y'.  Hit  <cont> 

2080  PRINT  "(To  the  left  along  the  flume) 

2090  PAUSE 

2100  CALL  Read_pot<"Y",V> 

2110  Y__m  i  n*  Y_z e  ro  + Y_s  iope*U 

2120  PRINT  "4.  Move  the  milling  machine  to  the 

2130  PRINT  "  maximum  value  of  ' y ' .  Hit  <cont> 

2140  PAUSE 

2150  CALL  Read_pot("Y",V) 

2160  Y _max  * Y_ze r o ♦  Y_s  iope«(7 

2170  PRINT  "5.  Move  the  milling  machine  to  the 
2180  PRINT  "  minimum  value  of  z  .  Hit  <cont> 

2190  PRINT  "(Downwards)" 

2200  PAUSE 

2210  CALL  Read_po t < "Z" , V) 

2220  Z _m i n* Z_ ze ro ♦ Z_s 1 ope»V 

2230  PRINT  "6.  Move  the  millinq  machine  to  the 
2240  PRINT  "  maxi mum  va 1 ue  of  z  .  Hit  <cont> 

2250  PAUSE 

2260  CALL  Read_pot("Z”,V> 

2270  Z_max ■ Z_ze ro » Z  slope'V 

2280  ! 

2290  1  E.  Write  out  the  file  ”Motor_coef " . 

2300  ’ 

2310  ON  ERROR  GOTO  Purge_file 
2320  Reenter:  CREATE  8DAT  Filenames,' 

2330  ASSIGN  «File  TO  Filenames 
2340  ' 

2350  OUTPUT  ?F i le : X_zero , X_s 1  ope 
2360  OUTPUT  WE i le ; Y_zero . Y _s 1  ope 

2370  OUTPUT  »Fi  le  :  Z_zero,Z_slope 

2380  ’ 

2390  OUTPUT  «F i le ; X_m i n , X_max 
2400  OUTPUT  WE  i  1  e ; Y'  min, Y_nax 

24  10  IJIJTPUT  «F  i  1  e  ;  Z_m  i  n  ,  Z_max 


2420 
243(3 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2520 
2530 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2630 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
27  70 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 


ASSIGN  «File  TQ  . 

SUBEXIT 

Purge_f lie:  PURGE  Filenames 
GOTO  Reenter 

i 

SUBEND 


SUB  Read_po t < D i rec t i onS , Va! ue> 

i 

i 

!  Read  one  potentiometer  and  return  a  volt 
'  age . 

• 

!  R3  --  Pot  X 

<  R4  --  Pot  r 

!  R5  --  Pot  Z 


IF  D i rec t i n$-''X"  THEN  Relay*3 
IF  Di rect ionS-"Y“  THEN  Relay-4 
IF  Di  rect  ion$-"Z"  THEN  Relay-5 


OUTPUT  723 : "OP . 1 ,0T” 

OUTPUT  723:"0B, 1 , " ;Re lay . 1 T“ 

t 

OUTPUT  723 : "IP ,  3T” 

ENTER  72301 -.Value 


SUBEND 


SUB  Clear_screen 
!  Clear  the  CRT, 

t 

OUTPUT  2  USING  •V,B":255. 75 

gclear 

SUBEND 


SUB  Mo  to r ( D i rec t lonS , Ro tat ion$ > 

i 

\  Turn  on  the  motor  in  the  requested  direc- 
1  tion  with  the  request  ion  rota- 

1  t ion  <CH.  CCH) . 


7  4 


3020  ’ 

’020  D i rS-D i rect ioni 

3040  RotS*Ao tat ioni 

2050  IF  Q  i  i"5“"X'‘  AND  RotS-"CU"  THEN  Lbit-2 
30S0  IF  Dir5-"X"  AND  RotS-,,CCW‘  THEN  Lbit-1 

3070  IF  D  i  rS-‘*Y"  AND  RotS-"CU"  THEN  Lbit-4 

3080  IF  D i rS-"Y"  AND  Rot$-"CCW  THEN  Lbit-3 

3090  IF  -D  i  r5-"Z"  AND  RotS-"CCW  THEN  Lbit-5 

3100  IF  D  i  r$»"Z"  AND  RotS-"CW  THEN  Lbit-6 

3110  ! 

3120  1 

3130  OUTPUT  723 : "OP , 7 , ” ;2*Lb i t ; ”T" 

3140  SUBEND 
3150  ! 

3160  ! 

3 1  70  * 

3  ‘  80  1 

3190  ! 

3200  ' 

3210  SUB  Motor_stop 

3220  ' 

3230  1  Stop  all  motors! 

3240  * 

3250  OUTPUT  723 : "OP , 7 , 0T" 

3260  SUBEND 
3270  ! 

3280  ! 

3290  ' 

3300  ' 

3310  ' 

3320  ! 

3330  SUB  Re t r ieve_coef <Coef ( *> , F i lenameS) 

3340  1 

3350  OPTION  BASE  t 

3360  ! 

3370  !  Retrieve  the  potentiometer  calibration 

3380  !  coefficients  from  a  disk  file  called 

3390  !  "Ho to r_coef " .  Place  these  in  an  array. 

3400  ! 

3410  ASSIGN  SFile  TO  Filenames 

3420  ' 

3430  FOR  1-1  TO  12  STEP  2 
3440  ENTER  *File;Coef(I> .Coef <I*1> 

3450  NEXT  I 

3460  1 

3470  ASSIGN  SFile  TO  » 

3480  SUBEND 
3490  ! 

3500  ' 

3510  ! 

3520  ! 

3530  ' 

3540  ! 

3550  SUB  Label _po i n t < X , Y , Z .Symbo IS) 

3r60  ' 

3570  !  Label  a  point  on  the  Flume  diagram 

3560  !  using  the  symbol  specif  ted. 

3590  ' 

3600  PRINT  TA8XY ( T  .13) . "<  X , Y , Z)  inches  - 
3610  PRINT  USING  "ODD. ODD" :X.Y.Z 
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3620 

t 

3630 

» 

3640 

IF  SVmbo 1 S-"X"  THEN 

3650 

» 

3660 

MOVE  14*Y.24+Z 

3670 

INDUE  -.2. -.2 

3680 

IDRAW  0.  A 

2690 

IDRAW  ,4,0 

3700 

IDRAW  0  ,  -  ,  4 

3710 

IDRAU  -.4,0 

3720 

i 

3730 

MOVE  14+Y,11*X 

37  A0 

I  WIVE  -  .2, -.2 

3750 

IDRAW  0  ,  .  4 

3760 

IDRAW  .4,0 

3770 

IDRAW  0  .  -  ,  4 

3780 

IDRAW  -.4.0 

3790 

END  IF 

3800 

i 

3810 

IF  Symbols-"*"  THEN 

3820 

t 

3830 

MOVE  1 4  +  Y  ,  1 1 *X 

38A0 

IMOVE  0 , .2 

3850 

IDRAW  -.2, -.4 

3860 

IDRAW  .4.0 

3870 

IDRAW  -.2,  .4 

i 

3880 

3890 

MOVE  1 4*Y  ,24*Z 

3900 

1MDVE  0.  .2 

3910 

IDRAW  -.2, -.4 

3920 

IDRAW  .4.0 

3930 

IDRAW  -.2,-4 

39A0 

END  IF 

3950 

SUBEND 

3960 

3970 

3980 

3990 

A000 

A010 

A020 

SUB  Move_ldv_to< X , Y , Z , Length 

DEG 


OUTPUT  723  ;  "CC  ,  1  T*‘ 

OUTPUT  723  ;  "5F , 3 . 3 . 3 , 1 ,25,12T" 

PRINTER  IS  1 

ON  KEY  0  LABEL  '•ABORT"  CALL  Stop_ail 
OPTION  BASE  ' 

Tolerance- . 00b 

Move  the  probe  to  the  position  indicated 
in  (inches)  relative  to  the  nottle  tip. 

Dimensions  are  in  inches  and  degrees 

A. I  Load  in  the  calibration  coeff icients. 

X_»ero*CoeT ( I ) 

X  s  1  one  *  L.je  •  (  2  ) 


« 


< _l  e  ro'i'oe1'  <  3 ) 

Y  _s  t  ope-Coe  *  <  4  * 

Z_:ero-Coef  <  5  > 

Z_s lope -Coef  <  b  > 

X_m i n-Coet  (  7 ) 

X_max -Coe f ( 8  > 

Y_n i n -Coe  f  <  9  > 

Y _«i» -Coef  (  1  0 ) 

Z_m i n-Coe  f <  t ! ) 

Z_max -Coef (12) 

i 

1  A. 2  Move  the  prooe 

i 

CALL  P  robe_"io«  ( Ang  1  e_p  robe ) 

i 

!  3.  Does  (X.Y.Z)  lie  within  the  per- 

1  mitted  boundaries7 

i 

IF  X>X_max  OR  X<X_m  in  OR  Y>Y_max  OR  Y<Y_min  OR  Z>Z_max  OR  Z'Z_mn  ^EKi 

4400  BEEP  1  700  , .5 

4410  BEEP  2000  ,  .5 

4420  BEEP  1  "00  ,  .5 

4430  BEEP  2000,  .5 

4440  PRINT  "Desired  point  is  out  of  range!" 

4450  SIJBEXIT 

44G0  END  IF 

4470  ! 

4480  !  C.  Find  out  where  the  probe  is  now,  draw 

4430  !  the  flume  on  the  CRT,  and  label  the 

4500  !  desired  position. 

4510  • 

4520  !  1.  Sound  warning,  movement  immenent! 

4530  ! 

4540  CALL  Clear_screen 

4550  PRINT  TABXYt 1 , 12) ."MOVEMENT  OF  MILLING  MACHINE  IMMINENT!"" 

4580  CQR  [-1  70  4 

4570  3EEP  1200,.  1 

4580  BEEP  1  700  ,  .  1 

4590  BEEP  2200  . .  1 

4800  BEEP  2700,  .1 

4610  NEXT  I 

4620  CALL  Clear_screen 

4630  ! 

4640  ! 

4650  OUTPUT  723 ;  "CC , 1 T"  ! CLEAR  RELAY  CARD. 

4660  CALL  Pos l t ion < "X” , X_ac t ual , Val u_5hat t , Lenqt h_arm .Leng th _p robe . Ang i e_arm . 
Ang 1 e_p robe , Coef <•> ) 

4670  CALL  Pos  i  t  lonCY"  ,  Y_actual  ,  Vaiu_shaf  t ,  Leng  th_arm  ,  Leng  th  _p  rooe  ,  Ang  1  e_am  . 
Angle_probe,Coef  <  * ) > 

4630  CALL  Pos i t i on< "Z", Z_ac t ual , Va 1 u_shaf t , Leng t h_arm ■ Leng th _p robe . Ang i e _arn 
Angle_proDe,Coef  <•>  > 

4690  ! 

4700  ’ 

4710  CALL  Draw_flume 

4720  CALL  Labe  1  _po  i  n  t  <  X_ac  tua  1  ,  Y_ac  t ual  ,  Z_ac  tua l  .” 4-"  > 

4730  CALL  Labe ! _po i n t < X . Y , Z . "X" > 

4740  ' 

4750  '  D.  Move  each  motor  to  brmq  the  error 

4760  '  between  actual  and  desired  position 


4210 
4220 
4230 
424  0 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4  370 
4380 
4330 
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4/70  ' 

4  780  ' 

4780 

4800  X  node 


i n  to  tolerance  . 


X  _o  1  d  *  X  _ac  tual 
ode:  Xe r r o r - X - x _ac t ual 

IF  A8S < Xe r ro r > > To lerance  THEN 
IF  Xerror>0  THEN  Rot5-"CCW" 

IF  Xerror<0  THEN  Rot$-"CW" 

CALi_  Mo  tor  <  “X"  ,  Ro  t$  ) 

CALL  Pos  1 1 1 on( "X” , X_ac  t ual . 9a  1 u_shaf t,Length_arm,Length_probe,Angie_ar 


Angle_prooe,Coet  <  • 
4860  ( 

a  i  > 

4870 

4880  GO' 

4890  END 

4900  ' 

4910  ' 

4920  Y_ol< 

493b  Y  node:  Yerri 


CALL  Plot _pat h (  X_o  Id  ,  Y_ac  t ual  ,  Z_ac  t ual  ,  X_ac  t  ual  ,  Y_ac  t ual  ,  7 

X_old-X_actual 
GOTO  X_node 
END  IF 


Y_o ld*Y_actual 
ode:  Ye r ror * Y - Y_ac tual 

IF  A8S ( Ye r ror > > To  1 e ranee  THEN 
IF  Yerror>0  THEN  Ro t S-  "CCW" 

IF  Yerror<0  THEN  RotS^'CW" 

CALL  HotorC'Y" ,Rot$) 

CALL  Pos  1 1  ion< “Y" , Y_ac tual , 9aiu_shaf  t , Lengt h_arro , Leng th _probe , Ang  :  e  _af 


Angie  _p robe , Coe f <*> > 


4990  CALL  Plot_path<X  actual,  Y  o Id , Z_ac tual , X_ac tual , Y_ac tual . /  jc 

al  > 

5000  Y  old*Y_actual 

5010  GOTO  Y_node 

5020  END  IF 

5030  ’ 

5040  ' 

5050  Z_old-Z_actuai 

5060  Z_node  :  Ze r ro r - Z - Z_ac t ual 
5070  IF  ABS< Zerror >> Tolerance  THEN 

5080  IF  Zerror>0  THEN  RotS-"CW" 

5090  IF  ZerrortO  THEN  RotS-"CCU" 

5100  CALL  Motor < "Z" .Ro t$ ) 

5110  COLL  PositionC'Z", Z_ac t ual . 9a 1 u_shaf t , Leng t h_arm ,Length_probe.Angle_ur 
Anqle_probe.Coef ( *  > > 

5120  CALL  P lo t_path ( X_ac tual , Y_ac tual , Z_o 1 d , X_ac tual , Y_ac tual . 2  ac 

a  1  > 

5130  Z_o 1 d“ Z_ac t ual 

5'40  GOTO  Z_node 

5150  END  IF 


5190 

CALL  Motor  stop 

5 ’90  ' 

5200  • 

5210 

FHR  1-1 

TO  4 

5220 

SEE° 

2400  , 

5230 

3EEP 

4800, 

.2 

524  0 

NEX’  I 

5250  ' 

5310  1 
5320  ! 

5330  SUB  Pas t t i on< D i r ec t i on5 . Val ue , Vai u_shaf t , Leng t h_ar m , Leng t h _p rooe . Png  1 
Angie_probe,Coef(*> ) 

5340  ! 

5350  OPTION  BASE  ! 

53S0  DEO 
5370  ! 

5380  1  Return  the  position  (inches)  for  the 

5390  '  aporooriate  direction  relative  to  the 

5400  1  nozzle  tip. 

5410  1 

5420  CALL  Read_pot(DirectionS,Vol tage) 

5430  ! 

5440  X_zero-Coef ( 1 ) 

5450  X_s 1 ope-Coe f (2) 

5480  Y_ze ro*Coef ( 3 ) 

5470  Y_s 1 ope^Coef ( 4 > 

5480  Z_ze ro *Coef < 5 ) 

5490  Z  s  lope*Coef (6 > 

5500  ! 

5510  ! 

5520  D i rS-D i rec t i on$ 

5530  ! 

5540  IF  DirS-"X"  THEN 

5550  Valu_shaf  t -X_zero+X_sl ope* Vo 1 tage 

5560  Value-Valu_shaf  t -Leng  th_arm*CQS ( Ang le_arm ) *TAN(  45-Ang le_arm/2 . 0  > 

5570  END  IF 

5580  ! 

5590  IF  Di r$*"Y"  THEN 

5600  Val  u_shaf  t*Y_zero+Y_slope*VoI  tage 

5610  Da I ue« Val u_shaf t+Length_probe*<l  .  -  COS(Anqle_p  robe  > >  +Leng  th_arn»' 

Ang le_arm) 

5620  END  IF 

5630  1 

5640  IF  D i rS*"Z"  THEN 

5650  Val u_shaf  t=‘Z_zero  +  Z_sIope“Vol  tage 

5660  Val  lie* Val u_shaf t-Length_probe*SIN(Angl e_probe ) 

5670  END  IF 

5680  1 

5690  SUBEND 
5700  1 

5710  ! 

5720  ' 

5730  ' 

5740  1 

5750  ! 

5760  SUB  S top_a 1 1 

5770  '  STOP  ALL  MOTORS  AND  QUIT 

5780  ! 

5790  GCLEAR 

5800  CALL  Motor_stop 

5810  PRINT  "MOTOR  CONTROL  ABORTED! 1'" 

5820  PRINT  "(HIT  <C0NT>  TO  CONTINUE)" 

5830  PAUSE 

5840  SUBEND 


5590  ' 

5900  ' 

5910  SUB  Plot_path<X1 , Y! ,Z1 .X2.Y2.Z2) 

5920  ' 

5930  !  Plot  the  path  of  the  probe  on  the  flume 

5940  1  diagram  as  the  motors  move  the  bed. 

5950  ? 

5960  !  Lower  plot  followed  by  upper  plot. 

5970  ! 

5980  PRINT  TABXY (l,t3>.”(X,Y,Z>  inches  - 
5990  PRINT  USING  "DO . ODD" ; X2 , Y2 , Z2 
6000  MOVE  14  +  Y1  .24  +  2! 

6010  DRAW  14+Y2.24+Z2 

6020  ! 

6030  MOVE  14+Yt , 1 1 +X! 

6040  DRAW  14+Y2.11+X2 

6050  ! 

6060  SUBEND 


I 


4.  T  SUBS 


a 


i  ib  j 

*•  o  i 

:o  ' 

-*0  SUB  T_coup 1 e<  Temperature,  Choices,  Scales,  No_re  ad  mgs,  S  t cev  > 
SO 
60 
7  0 
SO 
90 
100 
1  1  0 
1 20 
130 
1  40 
150 
160 
1  70 
130 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
230 
300 
310 
320 
230 
340 
350 


260 

OUTPUT  723:"CC,!T" 

370 

OUTPUT  723 ; "SF . 3 , 3 , 3 , 1 .25.12T" 

380 

♦ 

390 

OUTPUT  723: "OP, T ,0T" 

400 

» 

410 

f 

2.  Close  the  chosen  relays. 

420 

» 

430 

» 

a.  Ambient  T  --  Type  T,  relays 

6.9 

440 

1 

O.  Nozzle  T  --  Type  T,  relays 

7,8 

450 

t 

c.  Probe  T  --  Type  E,  relay 

9. 

460 

t 

470 

TF  ScaleS-"H"  THEN 

4  80 

'  Plot  a  histogram  using  Felsi 

is 

490 

HistoqramS“"YES" 

500 

ScaleS*“C" 

510 

ELSE 

520 

H i s  togramS-"N0" 

530 

END  IF 

540 

i 

550 

» 

560 

IF  Cho i ce$-"AMB IEN T"  THEN 

5  70 

T y  pe  5 ■ " T" 

5.-:  0 

OUTPUT  723  :  "OB , 1 ,6. : .3. 1 T” 

590 

END  1” 

-..in 

» 

SUBPROGRAM  T_C0UPLE 

by  Bill  Culbreth 

19  April  1984 

PURPOSE:  This  program  is  designed  to 

read  type  T  or  E  thermocouples  and 
return  the  actual  temperature  to 
the  calling  routine. 

Temperature  --  Temperature  from  the 

t he rmocoup 1 e  in  degrees  F  or  C  from 
the  t he rmocoup le  identified  by 
"Cho  iceS" . 

Choices  --  Thermocouple  choice,  current¬ 
ly:  "AMBIENT",  "NOZZLE”,  or.  "PROBE". 

ScaieS  --  "F"  for  Fahrenheit  or  "C"  for 
Celsius,  or  "H"  for  histogram  in  "C" . 

No_readinqs  --  Hou  many  readings  of  the 
same  thermocouple  should  tne  routine 
take? 

Stdev  --  The  standard  deviation  of  the 
temperature  for  the  indicated  numoer 
of  readings  in  the  units  given  by 
"Scales” . 


1.  Open  all  relays  and  initialize  the 
A/D  converter. 


■ 
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I 


B.  :  . 


END 

» 

IF  2ho iceS’‘”DRCBE"  'HEN 
Ty oe$ - "E" 

GUTP'.i 7  ’23 :  "08 , 1  .3,  l  T 
END  IF 


2.  Taire  an  A/D  conversion  and  convert 
into  temperature. 

Sum*0 
3  um 1  -  0 

HA  IT  1 

IF  H j s tog  rami-" VES"  THEN  GOSUB  Set_up_Nisto 

t 


FOR  I-l  TO  No_readinqs 
OUTPUT  723  ;  ■*  IP  ,  JT" 

ENTER  72201  ;  A 

IF  'ype'S*'1 T"  THEN 
A-A/1Q00 
END  t.p 

IF  TypeS* "E"  THEN 

A  *  A  / ' r!  0  0 
CND  IF 

P  I N  T  "V<nW>  •  "  ;  A 
GOSUB  Convert  t 

BEE:’  'A»1Q0,  .01 
3um*Sum* A 
Sum  1 'Sum  1  ♦  A*A 

t 

IF  H i s t og ramS* "TES"  rHEN  GOSUB  p!ot _?oint 
NE  <  ’  I 

remperature*Sum/No_readings 


IP  No_read i ngs* 1  'HEN 

j tdev  *  0 

else 

S tde v-SQR ( ABS < < Sum  1 -No_readinqs*'emperature  2>/<No_readmg 

end  IP 

A .  Open  all  re  1  ays  . 

TlJrpllT  723  ;  "OP  .  1  ,  0T" 

IJBEXIT 


nv**  r  t  _  t  : 


suDroutine  conver • s  ( nv ) 


1210  '  erature  . 

1220  ' 

1230  ' 

1240  IF  TypeS-"T"  THEM 

1250  A-2.5bS129  7E*1*A-S.  1 954869E-  1  *A*A*2 . 2 1 8 1  644E-2*A '  3-3 . 55009E-4-A  *  4 

1250  END  IF 
1270  ! 

1280  r 

1290  IF  T ypeS-"E"  THEN 

1300  A-1  .  7022525E  *-1  *A-2.2097240E-1*A«A +5. 480931  4E-3»A~3~5. 76G9892E-  5*8  ’  4 

1310  END  IF 
1320  ! 

1330  !  Fix  the  scale. 

1340  * 

1350  IF  ScaleS-"F"  THEN  A-1.8*A*32 
1350  RETURN 
1370  ! 

1380  1 

1390  ! 

1400  ' 

1410  * 

1  420  ' 

1430  Set_up_h isto :  ! 

1 440  t 

1450  1  Set  up  a  histogram  of  temperature 

1460  1  versus  numDer  of  counts. 

1470  ! 

1480  !  1.  Zero  out  the  Height(*>  array. 

1490  ' 

1500  0 rn  He tght< 203) 

1510  ' 

1520  FOR  I-i  TO  202 

1530  He  l  gh  t  <  I )  -0 

1540  NEXT  I 

1550 

1560  ' 

1570  DUMP  DEVICE  IS  701 

1530  GIN  IT 

1590  OUTPUT  2  USING  "*  ,  B”  :  255 . 75 

1600  GRAPHICS  ON 

1610  FRAME 

1520  WINDOW  -100,100,-10.100 

1630  MOVE  -65.82 

1640  OSIZE  7 

1650  LABEL  "TEMPERATURE  HISTOGRAM" 

1660  AXES  25,10.0,0.4,5.3 

1670  PEN  -1 

1660  MOVE  0,-10 

1630  DRAM  0,0 

1700  MDVE  0,00 

1710  DRAW  0.  '00 

1 720  PEN  ’ 

1  730 
1  740 


lake  10  temperature  readings  to  get 


1810  ENTER  72301 :A 

1 320  NEXT  [ 

1330  ' 

1 840  FOR  1-1  TO  1  00 

1850  OUTPUT  723 ; "IP , 3T" 

I860  ENTER  72301 :A 

1370  IE  Type’S-" T"  THEN 

1880  A-A/1000 

1830  END  IF 

1900  IF  TypeS-"E’*  THEN 

1910  A-A/1000 

1920  END  IF 

1930  GGSU8  Convert_t 

1940  Sum-Sump-A 

1950  Sum  1 -Sum  1  *A*A 

1960  NEXT  I 

1970  T_mean-Sum/1 00 

1980  Sd-SQR<A8S(Sum1-100*T_mean  2>/99> 

1990  ' 

2000  !  Change  the  window  to  extent  from 


2010 

'  -4-Sd 

to  *4»Sd. 

2020 

i 

2030 

MOVE  -12.-7 

2040 

CSIZE  4 

2050 

LABEL  USING 

"ODD . DD“ ; T_mean 

2060 

MOVE  15,-7 

2070 

LABEL  "C" 

2080 

i 

2090 

MOVE  -87,-7 

21  00 

LABEL  USING 

"DDD . DD" ; T_mean-3*Sd 

21  10 

MOVE  63,-7 

2120 

LABEL  USING 

"DDD  .  DD-"  ;  T_mean  +  3*Sd 

2130 

i 

2140 

i 

2150 

MOVE  5.47 

2160 

LABEL  "50" 

2170 

* 

2180 

i 

2190 

1  3.  Calculate  the  window  temper 

2200 

!  ature 

interval . 

2210 

♦ 

2220  I n te rva l -4*Sd/ 1 00 

2230  ' 

2240  T_mi n* T_mean-4«Sd 

2250  T _max  *  I _mean»A*Sd 

2260  ! 

2270  WINDOW  T_min,T  max, -10, 100 

2280  RETURN 
2230  * 

2300  ' 

2310  ' 

2320  ' 

2330  ' 

2340  ' 

2350  PI o t_po i n t :  ! 

2360 

2370  '  Plot  each  temperature  point  as 

2  S8D  '  received  on  the  histogram. 

2330  ' 

JO0  -req-((A-T_minl/(T_max-^_min))»3000*-1000 
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241 0  BEEP  F  req  .  .01 

2420  J-0 

2430  ' 

2440  FDR  T»T_min  TO  T_max  STEP  Interval 

2450  J-J  +  1 

2460  IF  A<-T  THEN 

24/0  He i qh t ( J ) “He i gh t < J > + 1 

2480  GOTO  Continue 

2480  END  IF 

2500  NEXT  T 

2510  ! 

2520  1  Draw  the  point. 

2530  ' 

2540  Continue:  1 

2550  1 

2560  ' 

2570  MOVE  A  ,  He  lgh  t  <  J  > 

2560  CSIZE  4 

2590  LABEL  'V 

2600  ’ 

2610  ?  Label  the  number  of  points  and 

2620  1  the  temperature. 

2630  1 

2640  PRINT  TABXY<2 ,3) ;“T ( C > - ”  ;A 

2650  PRINT  TABXY ( 2 ,5) : “Sample  *";I 

2660  1 

2670  RETURN 
2680  ' 

2690  f 
2700  ' 

2710  ’ 

2720  ! 

2730  1 

2740  SUBENO 
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i  • 


5  .  T  CAL 


■o  •  t_:al 

:o  ’ 

50  1  rhis  proqram  records  the  output  t  rom  the 

40  lamoient  ,  proDe  and  no;  ;  le  tnsmoitoupies 
50  'for  calibration  purposes.  'he  data  is 

SO  'printed  as  uie  l  i  as  transmitted  to  a  data 

,70  'file. 

00  1 

30  (Data  is  recorded  in  the  following  order: 

’00  '  a.  Ambient  temperature 

l(Q  (  b.  Probe  temperature 

120  1  c.  Nozzle  temperature 

130  ( 

140  'Place  the  ambient,  probe  and  nozzle 

150  (thermocouples  into  a  bucket  of  warm  water, 

160  (execute  this  program,  then  add  ice. 

170  ’ 

130  * 

130  LOADSUB  ALL  FROM  ••T_SUBS•• 

200  ( 

210  !  Identify  the  BOAT  file 

220  ! 

230  INPUT  "NAME  OF  FILE  WHERE  DATA  IS  TQ  BE  STORED?" . F i lename 1 S 
240  ! 

250  !  Identify  the  number  of  data  points  desired 

260  ’ 

270  INPUT  ''NUMBER  OF  DATA  POINTS  DESIRED?", N Items 
280  ( 

230  '  Identify  the  number  of  samples  per 

300  ?  thermocouple  per  data  point 
310  ( 

320  INPUT  "NUMBER  OF  SAMPLES  PER  DATA  PO INT?" , Samp le 

330  ' 

340  !  Initialize  a  counter 
350  ( 

360  1-1 

370  ' 

330  (  Create  a  data  file 

3  30  ' 

400  Recor ds* < N i t ems*8*3/256 ) *2 
410  CREATE  BOAT  F  i  i ename 1 S , Reco rds 
420  ASSIGN  SFilel  TO  Filename  IS 
430  ' 

440  !  Taxe  data  and  print  results 

450  ' 

460  CALL  T  coup  1 e < T_amb i ent , "AMB IEMT" , "C" , Samp  1 e , Sd ) 

470  BEE?  T *50  .  .05 

430  PRINT  "TAMB , SD-": T_ambient ,Sd 

490  CALL  T_coupie< T, "PROBE”, "C",3ample,St_dev> 

500  BEEP  T *50  .  . 05 

510  PRINT  "TPR0BE.SD-": T.St_dev 

520  CALL  T_coup i e f T_no 1 1 ie . "NOZZLE" , "C” , Samp  1 e , Sd ) 

530  BEEP  T-50..05 

540  PRINT  "TNUZL.SD-";T  nozzle. Sd 

550  ! 

560  !  Send  data  to  the  data  file 

570  ' 

580  OUTPUT  -F  i  le 1  :  T  ambtent.T.T  nozzle 

590  ' 

p00  !  Test  to  see  it  finished 
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bio  ' 

620  IF  I*Ni terns  THEM 
630  GOTO  760 

640  ELSE 
650  I  *  I  ♦ 1 

660  bfl  i  0  460 

670  END  IF 
680  ' 

630  !  Close  the  data  file 

700  ! 

710  OUTPUT  iFilet : -  1 00 
720  ASSIGN  if i lei  T0  * 

730  ! 

740  1  Alert  the  user.. .the  job  is  completed 

750  ' 

760  PRINT  "ALL  DONE" 

770  BEEP 
780  BEEP 
730  END 


PROBE  CAL 


:  0 

_0 
20 
40 
50 
60 
70 
80 
30 
!00 
1  1  0 
120 
130 
140 
150 
160 
170 
180 
130 
200 
210 
220 
230 
240 
250 
26U 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
330 
390 
400 
410 
420 
a30 
440 
450 
460 
470 
430 
430 
500 
510 
520 
530 
540 
550 
560 
5  70 
5.30 
530 
600 


--RuBE  _lAl 

'his  proqram  will  provide  the  data 
necessary  to  :ai  iDrate  prooe  oef.ec'.on 
ii  used  in  the  following  manner: 

a.  Switch  OFF  probe  actuation  power 
at  the  probe  Dase . 

b.  Swing  the  probe  arm  out  of  the 
tan*  such  that  the  proDe  is 
accessab le . 

c.  Attach  the  calibration  panel  to  the 
probe  assembly,  taxing  care  not  to 
damage  the  glass  proDe. 

d.  Run  the  proqram... it  will  ask  for 
the  desired  position  in  mV.  The 
following  guidelines  apply: 

l.  If  it  is  desired  to  lower 
the  probe,  type  the  extreme 
value  9000  (anything  >4600 
will  work,  but  9000  is  a 
quick  and  easy  numDer  to 
enter) . 

ii.  If  it  is  desired  to  raise 

the  Probe,  enter  the  extreme 
value  90  (anything  less  than 
940  will  wo  rx  >  . 

e.  The  program  will  next  ask  which  bit 
is  selected  to  be  ’‘high".  The 
following  guidelines  apply: 

i.  If  it  is  desired  to  lower  the 
probe .  enter  7 . 

ii.  If  it  is  desired  to  raise  the 
probe,  enter  3. 

f.  Switch  ON  probe  activation  power  and 
when  the  probe  reaches  a  desired 
degree  of  deflection,  switch  the 
activation  power  OFF  and  record  the 
mV  value  printed  on  the  screen. 

g.  Repeat  steps  <d>  through  (f)  until 
sufficient  data  is  collected.  Load 
this  data  into  the  program  "P0LYFII" 
and  request  a  second  order  fit 

(let  X»anqle  and  Y-mV>. 

Enter  the  coefficients  derived  by 
this  urogram  into  the  appropriate 
location  m  tne  proqram  "DRDBE_SI IBS” 


NOTE: 


The  relays  will  be  activated  on 
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MOTOR  CAL 


10  '  '10T0R_CAL 

JO  ' 

50  !  ^'JRPQSE :  Calibrate  the  potentiometers 

AO  '  used  with  the  motors  on  the  milling 

50  '  macnine. 

50  ' 

70  ' 

30  OP T ION  BASE  ! 

90  DIM  Coef  <  1  2 > 

100  ' 

110  OUTPUT  723 ; "SF .3.3,3.! .25 . 1 2T" 

120  LOAOSUB  ALL  PROM  "MTR_SUBS" 

130  L0ADSUB  ALL  PROM  "°ROBE_SUBS" 

1A0  P 1 1 "mo  tor _coef " 

150  ' 

160  CALL  Cal i b rate < F i le$ > 

1 70  ’ 

180  CALL  Re t r t eve_coef < Coet < * > , F i 1 eS ) 

190  BEEP  2A0O  .  1 
200  1 

210  1 

220  PRINT  "What  (x.y.z)  position  do  you  wish  to" 

230  PRINT  ”  move  toR  (inches  relative  to  noz-" 

2 A0  PRINT  "  z 1 e > “ 

250  INPUT  X.Y.2 

260  ! 

270  CALL  Move_ldv_to< X , Y , Z ,Length_arm .Length _probe . Ang le„arm , Probe_ang i e , Coe'  < 
• )  ) 

280  ' 

290  BEEP 

300  GOTO  290 

310  END 


90 


LOAD  XYZ 


'  i_0AD_<YZ 

t 

!  by  Bill  Culb re th 

> 

'  30  April  !984 


!  PURPOSE:  This  program  will  allow  the 
!  user  to  enter  desired  <x.y,z>  posi- 

!  tions  of  the  milling  machine  relative 

!  to  the  nozzle,  The  values  will  be 

!  stored  on  disk  to  De  utilized  later 

*  by  M A I N _ T . 

i 

DIM  X<  500 )  ,Y<500> ,Z<500> 

i 

?  !.  Input  the  file  name. 

t 

GOSUB  Clear  screen 


PRINT  "Input  the  file  name." 

PRINT  ”{  I  suggest  'RUNXX'  where  'XX” 
PRINT  "  is  the  run  number. 

PRINT 

INPU T  F i lenameS 


Begin  inputting  data. 


GOSUB  Clear_screen 

i 

'RINT  "Do  you  wish  to  append  a  previous  data  file? 
INPUT  Pnswer'i 

i 

IF  AnswerS-"YES"  THEN 

INPUT  "Previous  file  name7", Old  file£ 

ASSIGN  iFilel  TO  Old  file's 


1-1*1 

ENTER  «Fi  lei ;X( I) , Y( I) . Z< I) 

OUTPUT  701:"I,X.Y,Z*";I.X<I).Y(I),Z<I) 
IF  X<IX>-I00  THEN  GOTO  Loopl 
ASSIGN  iFilel  TO  - 
PURGE  0 1 d_f i 1 eS 
Conn  t  *  I  -  1 
ELSE 

Coun  t - n 
END  IF 


GOSUB  CIear_screen 

BEEP 

I 

PRINT  "1.  Input  the  desired  position  i 
PR  IN'  "  mcnes  as  X,  Y ,  ano  I." 

PRINT 

PRINT  "2.  'erminate  input  by  entering 


610 

PRINT 

••  n 

’ernmate  input  by  entering  ’-'NO"" 

620 

PRINT 

** 

for  X  , 

r  ,  ana  2." 

630 

PRINT 

"3 

If  you 

wish  to  enter  an  orientation  angle." 

640 

PRINT 

" 

enter 

-993, orientation, 0  . " 

550 

PRINT 

S60 

» 

670 

Begin:  1 

680 

Coun  t -Coun t  +  1 

690 

PRINT 

'Item  *" 

:Count 

700 

INPUT 

‘  <  X  ,  Y  ,  Z ) 

in  inches 7 “,X(Count> ,Y(Count> , Z ( Coun  t 

710 

OUT 

PUT 

701  :  "I , 

X.Y,Z*";Count.X<Count> , Y ( Coun  t ) , Z ( Coun 

720 

IF 

X ( Coun  t ) <  >- 

100  THEN  GOTO  Begin 

730 

1 

740 

End 

coun  t  »Coun 

t-t 

750 

» 

760 

All 

data  po 

ints  have  been  entered. 

770 

♦ 

730 

t 

a . 

Set  up 

new  sof  tk  eys . 

790 

» 

b  * 

Exp  Lam 

sof  tk  ey  s  . 

800 

t 

810 

GOSUB  Clea 

r  ..screen 

820 

i 

830 

PRINT  1 

SOF 

TKEY  LABELS:" 

840 

PRINT  ' 

0 

--  WRITE  data  to  disk." 

850 

PRINT  • 

2 

--  EDIT 

out  bad  data." 

860 

PRINT  ’ 

4 

--  HARD 

copy  the  data  on  printer." 

370 

PRINT  • 

6 

--  LIST 

data  on  the  CRT ." 

880 

PRINT  ' 

8 

--  STOP 

terminates  the  program." 

890  ' 

900  ON  KEY  0  LABEL  "WRITE"  GOSUB  Wnte.data 

910  ON  KEY  2  LABEL  "EDIT"  GOSUB  Edit_data 

920  ON  KEY  4  LABEL  "HARD"  GOSUB  Hard-copy 

930  ON  KEY  6  LABEL  "LIST"  GOSUB  List.data 

940  ON  KEY  3  LABEL  "STOP"  GOTO  Terminate 

9SQ  GOTO  900 

360  f 

970  ' 

380  ' 

990  ! 

1000  ' 

1010  ' 

1020  1 

1030  Clear_screen:  ! 

1  040  1 

1050  1  Clear  the  CRT  display. 

1  060  1 

1070  OUTPUT  2  USING  , B" : 255 . 75 

1080  GCLEAR 

1090  RETURN 
1100  ' 

1110  1 

1120  ' 

1130  1 

1140  ' 

I  150  ' 

HSO  Hard_copy  :  ' 

1170  ’ 

1190  '  Print  out  all  data  to  the  orinter. 

’190  ' 

'  200  PRINTER  IS  701 
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1 2 1 0  ' 

'220  PR  TNT  "  Count  X(  inches)  /<count>  Z<count)’ 

1230  PRINT  "  _  _  _  _ ' 

! 240  PRINT 

1250  FOR  1=1  TO  £nd_count 

'260  PRINT  I .  X(  I )  ,  Y<  I  >  ,  Z<  I  > 

1270  NEXT  I 

1280  1 

1290  PRINTER  IS  1 

1300  RETURN 

1310  1 

1320  1 

1330  1 

1340  1 

1350  • 

1360  ’ 

1370  List_data:  1 
1380  ' 

1390  !  List  data  to  the  CRT. 

1 400  ! 

1410  G0SU8  Clear_screen 

1420  ! 

1430  PRINT  "There  are" ;£nd_count ; "  data  points.” 

1440  PRINT 

1450  INPUT  "Uhich  point  do  I  start  unthR" .Start 

1460  INPUT  "which  point  do  I  end  u i t h7" , End_data 

1470  1 

1480  FOR  I-Start  TO  End_data 

1490  PRINT  " I , X , Y , Z * " : I , X ( I >  ,Y(I)  ,Z(I> 

1500  NEXT  I 

1510  RETURN 
1520  1 

1530  1 

1540  1 

1550  ' 

1560  1 

1570  1 

1580  Edit_data:  1 
1590 


1600 

?  Edit  data. 

1610 

t 

1620 

G  OSUB 

C Iear_sc  reen 

1630 

1 

1640 

INPUT 

"Which  data  point 

do  you 

wish  to  al ter7" , I 

1650 

» 

1660 

PRINT 

1670 

PR  TNT 

"For  point 

,  <X  ,Y  , 

Z)  where:" 

1680 

PRINT 

X<  I)  ,  "  ;  Y'<  I)  :  "  ,  " 

;Z<  D 

1690 

t 

1  700 

TNPIJT 

"Type  in  the  new 

va l ues : 

' ,  X ( I ) , Y ( I ) ,Z(I) 

1710 

RETURN 

1  720 

» 

1  730 

1 

1  740 

t 

1  750 

♦ 

1  760 

♦ 

1  770 

t 

1  780 

Ar  i  te_Uata : 

» 

1  790 

» 

'800 

'  dr 

i  te  data  out  to  a 

tile  on 

disk. 

9  3 


:  8  ’ 1  j  ' 

'  3 0  <  <  Ena_count  *  i  >  *  •  l  00 

"830  v  <  End  _coun  t  00 

'340  Z  <  E nd  cijun  t  ♦  1  *  ■  -  l  00 

V350  ’ 

'860  Max_da ta* 3* < End_coun t ♦ I > 

1370 

i860  File_site*rNT<Max_data*8/256>+1 

1880  ' 

1900  CREATE  9DAT  F  i  lenameS  .Fi  le_s  i  ze 

1910  ASSIGN  '«File  TO  Filenames 

1920  ' 

1930  FOR  I*'  TO  End_cou.n t  +  1 

1  940  OUTPUT  '*F  i  J  e  ;  X <  I )  ,  Y  <  I  > .  Z <  I  > 

1950  NEXT  I 

1960  1 

1970  G0SU8  Clear_screen 

1980  BEEP  2400.  .3 

1990  PRINT  "File  " ;F i lenameS ;  '*  has  been  stored! 

2000  ASSIGN  *F  i  le  TO  * 

2010  RETURN 
2020  ' 

2030  ' 

2040  ! 

2050  ! 

2060  ' 

2070  ! 

2080  Term i nate :  ' 

2090  G0SIJ8  Clear_screen 

2100  PRINT  "NORMAL  TERMINATION  OF  PROGRAM!" 

2110  END 


94 


SEND  DATA 


9  . 


■  0 
20 

30 
40 
50 
60 
70 
80 
'30 
1  00 
1  !  0 
120 
130 
1  40 
150 

160  Sc-9  ’  RS-232  IS  SELECT  CODE  9. 


1  70 

PRINTER  IS  1  ' 

PRINTER  IS  CRT. 

180 

Pr- 1  ' 

DEFAULT  PRINTER  IS  i.RT 

130 

P  r i n  t er _cho ice* 

701  ’  MY  PRINTER  IS  701 . 

200 

3 i ts*  7  1 

BITS  PER  CHARACTER 

210 

Duplex*0  ' 

full  duplex 

220 

8aud*300  1 

BAUD  RATE 

230 

Computer*'  1 

ASSUME  IBM  COMPUTER 

240 

i 

250 

OUTPUT  Pr  ;  " ( 300 

BAUD,  IBM  assumed." 

260 

OUTPUT  Pr: "  Load  the  binary  prop  ran  3EB  •  :r 

270 

OUTPUT  Pr;"  unless  you  have  BASIC  2.3' 

280 

OUTPUT  Pr;"  SET 

MODEM  ON  <FULl  DljP.t'o  " 

290 

OUTPUT  Pr;"  " 

300 

i 

210 

DIM  NameST 200 1 , Ho  f i le$( 30 ) , Aa< 1 500 > . NiunoSt 

320 

IN  fEGER  Iseno 

330 

» 

340 

CONTROL  Sc. 3; Baud 

350 

CONTROL  Sc,4:Bits-5+4  '  BITS/CHAR  1  *STCP  3 

360 

i 

370 

i 

280 

To_d isk  *0 

230 

0ataaumo*0 

400 

I  data- 1 

41  0 

1*1 

420 

J-l 

4  30 

<  *  1 

44(0 

L  *  1 

450 

t 

460 

ON  ERROR  GOTO  E 

r  ro  rs 

h  70 

ON  kEY  0  l ABEL 

"Line  Mode"  GOTO  Line_mode 

480 

ON  *EY  5  LABEL 

"Terminai"  GOTO  Terminal 

490 

ON  “'EY  6  L'-RE- 

"To  Or:"  GOTO  Pr  crt 

500 

ON  «'EY  7  CARE- 

"To  Pr  t"  GOTO  ?r_prt 

510 

ON  LEY  8  i_ABEL 

"DATA"  GOTO  Data_.'ump 

520 

1 

530 

t 

54  0 

Line  mode ;  ’ 

550 

OUTPUT  Pr; 

"'LINE  RECEPTION  MODE  1 " 

560 

3eq  in:  S  TA  ’!JS  Sc  . 

1  0  :  Y  •  CHECK  "OR  “  * : L BUFFER 

570* 

580 

5  50 

IF  3 1  T ( Y  ,  0  )  - 0  THEN  GOTO  Begin 

-.01) 

•  RECEIVE  ROUT 

INF. 

SEND_DATA 

To  VAX. I3M. IRS-90. 

HP-9826  TERMINAL  PROGRAM 
[REQUIRES  BINARY  ENHANCEMENT  PROGRAM 
"8E8" '  ] 

JUNE  30.  1982 
updated  1 /5/83 
updated  I / i 6/84 

3 ILL  CUL3RETH 
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HD-A152  224 


UNCLASSIFIED 


DESIGN  AND  CONSTRUCTION  OF  A  COMPUTER  CONTROLLED 
HIROTHERMOCOUPLE  PROBE  FOR  THE  STUDV  OF  BUOVANT  JETS 
(U)  NAVAL  POSTGRADUATE  SCHOOL  HONTEREV  CR 
R  J  HATOUSHEK  SEP  84  F/G  9/2 


NL 


p 


MICROCOPY  RE-SOLUIlON  TLST  CHARI 

V.'.riN.V  ' 


•••■N!  >.u-i  •  ; 


bi  o  r 

620  Recewe; 
630  B-A 


630  B-A 

bAO  OUTPUT 

650  IF  B-6: 

bbO  I F  B  - 1 

6/0  GOTO  8. 

680  * 

690  !  TRANSM 

700  ' 

710- 
720 
730 
7  AO 
750 
760 
770 
780 
790 
800 
810 
820 
830 
8A0 
850 
860 
870 
880 
890  ' 

900  • 

910  ! 

920  '  DATA  F 

930  1 

9  AO  1 
950  ' 

960  Oata_dump: 
970 
980 
990 
1  000 
1010 
1020 
1030 

I  0A0 
1050 
1  060 
1070 
1080 
1090 
1100 

1 i 1 0  GOTO  Beg  i  n 
1120  ' 

M  30  ' 

I I  AO  '  ERROR  H 
1150  ' 

1160  Errors:  G 


STATUS  Sc  . 6  ;  A 


OUTPUT  or  USING  .  A‘* :  CHR$<  B > 

IF  B-b3  AND  0atadump-1  THEN  GOTO  Data_dump 
IF  B - 1 3  AND  Computer-3  THEN  OUTPUT  Pr:CHRS(13> 
GOTO  Begin 

TRANSMIT  ROUTINE. 


IF  Duple*-0  THEN 

IF  NUM<KeyS><>255  THEN  OUTPUT  Pr  USING  "# .  A" ; Keyi 
IF  NUM<  Key  $ ) *255  THEN  OUTPUT  Pr;"  " 

END  IF 

IF  Computer-1  AND  NUM(Key$>-8  THEN  Key$-CHR$( 6A ) 

i 

!  the  previous  line  gives  an  * 

!  for  a  backspace  for  tne  I3M. 

t 

IF  Computer-5  ANO  NUH<Key$>-8  THEN  KeyS-CHR$< 1 271 
» 

*  THE  VAX/VMS  REQUIRES  A  DELETE 
!  SYMBOL  FOR  A  BACKSPACE. 

t 

IF  NUM(KeyS) -255  THEN  Key$-CHRS< 1 3) 

OUTPUT  Sc  USING  "•  .A";KeyS 
GOTO  Begin 


DATA  FILE  OUT  TO  THE  HOST  COMPUTER. 


IF  I_data-1  THEN  GOSUB  Open_file 

i 

IF  Datadump-0  THEN  GOTO  3egin 
IF  Computer-1  THEN  HAIT  .3 
'  wait  for  the  slow  I3M. 

BEEP  1 000 *RND- 1 500 ,  .05 
OUTPUT  P r ; ” A ( ” ; I_da ta; " ) -" ; 
OUTPUT  Pr ; Aal I _da t a ) 

GOSUB  Send_numbe r 
IF  Aa< I _da t a) - -200  THEN 
I_data- 1 
Datadump-0 

END  IF 

I_data-I  data»l 


ERROR  HANDLING  SUBROUTINE 


OFF  ERROR 
C  1  file  *  t  i  e  -  -200 
'  r I RS T ,  END  3F  FILE  ERROR. 
IF  FRRN-59  THEN 
Aa< I )--2Q0 


GOTO  2000  '  RETURN  AFTER  ERROR. 

END  IF 

i 

IF  ERRN059  THEN  OUTPUT  Pr;"<error  *"  ;  ERRN  :  "  generated  .  >" 

IF  ERRN-54  THEN  OUTPUT  Pr ;  f  FILE  <";Hp_f iie$;">  ALREADY  THERE!" 
IF  ERRN-54  THEN  GOTO  Created 

IF  ERRN-56  THEN  OUTPUT  Pr;"<FILE  <“ ;Hp_f i Ie$ ; "  IS  NOT  ON  OIS< . >' 
ASSIGN  jFile  TO  * 

GOTO  Lme_mooe 

t 


1210  GOTO  2000  '  RE' 

1220  END  IF 

1230  * 

1240  IF  ERRN059  THEN  i 

1250  IF  ERRN-54  THEN  Ol 

1260  IF  ERRN-54  THEN  G( 

1 2 TO  IF  ERRN-56  THEN  Ql 

1280  ASSIGN  sFile  TO  * 

1290  GOTO  Lme_mooe 
1300  1 
1310  ! 

1320  1  OUTPUT  TO  CRT. 

1330  ! 

1340  P r_c r t :  P r *  1 

1350  GOTO  Line_mode 

1  360  * 

1370  * 

1380  1  OUTPUT  TO  PRINTER. 

1390  1 

1400  Pr_prt:  Pr-Printer_choice 
1410  GOTO  Line  mode 


!  CHANGE  THE  TERMINAL  CHARACTERISTICS. 

i 

Term i nal :  ’ 


1470 

OUTPUT 

Pr;" 

1  . 

Baud  Rate  *";Baud 

1480 

OUTPUT 

Pr;" 

2. 

Bits/Char  *":Bits 

1490 

OUTPUT 

Pr:" 

3. 

Duplex  •“{Duplex 

1500 

OUTPUT 

Pr;" 

[1-full .0-half 1" 

1510 

OUTPUT 

Pr;" 

4. 

Computer  »“;Computer 

1520 

OUTPUT 

Pr;" 

1  IBM-1  .  VAX/UNIX-2.  " 

1530 

OUTPUT 

Pr;" 

TRS-80-3,  Cyber*4,  v. 

1540 

OUTPUT 

Pr;"  " 

1550 

INPUT 

'Change 

which 

one?" .Which 

IF  Which-1  THEN  INPUT  "To?" . Baud 
IF  Which-2  THEN  INPUT  "To?". Bits 
IF  Which-3  THEN  INPUT  " To?" . Dup lex 
IF  Which-4  THEN  INPUT  "To?" .Computer 
IF  Computer*'  THEN  Duplex-0 
IF  Computer*3  THEN  Duplex-0 
IF  Computer*3  THEN  Bits-8 
IF  Computer*5  THEN  Duplex*! 


1640  GOTO  Line_mode 
1650  1 
1660  ! 

1670  ! 

1 680  Open_f i 1 e :  * 

1690  ! 

1  700  1 

1710  ! 

1 720  Datadump* 1 

1 7P0  ' 


Open  a  file  to  read  data  from 

d  i  sk  . 


INPUT  "Is  this  LDV  data’  ( 1 -YES)"  ,Ldv$ 

IF  Ldv$  - " 1 "  THEN 

INPUT  "Experiment  #?" . E xper i men t$ 

ELSE 

OUTPUT  Pr:"Data  f lie  out  of  HP  to  host." 
INPUT  "Fi[p  njnel",Hp  •lies 
END  IF 


97 


I  I.  J  I II 


"510  1 

1P20  IF  LdvS-'T-  THEN 

1830  Hp_f  i  leS“ExpertmentSa''_RESULT" 

1840  END  IF 

1 950  • 

I860  !  Read  the  file  off  of  disk. 

18  70  • 

1880  ASSIGN  if  ile  TO  Hp_f1le$ 

1890  1*1 

1900  Check *0 

1910  BEEP 

1920  BEEP 

1930  OUTPUT  P  r  ;  “  ( Work  i  ng  on  file  <**  ;  Hp_f  i  le$  ;  ’*>  . } 

1 990  1 

1950* 

1960  ENTER  4File:Aa<I> 

1  970  Check  *Aa< I ) 

1980  1-1*1 

1990* 

2000  • 

2010  ASSIGN  SFi le  TO  • 

2020  Datadmnp* 1 

2030  RETURN 
2040  1 
2050  * 

2060  > 

2070  Send_number :  1 
2080 
2090 
2100 
2*10 
2120 
2130 
2140 
2150 
2160 
21  70 
2180 
2190 
2200 
2210 
2220 
2230 
2240 

2250  RETURN 
2260  ' 

2270  ' 

2280  ' 

2290  END 


'  SEND  A  NUMBER  ONE  CHARACTER  AT 
1  A  TIME  TO  THE  HOST  COMPUTER. 

t 

Numb$«VAl$<Aa< I_data>  > 
Length«L£N<Niunb$) 

i 

IF  <<LdvS- "!■•>  AND  ( I_data>  1  3>  >  THEN 
Posit*POS<Numb$."."> 

IF  (Posi  tOO)  THEN  Length*Pos i t+2 
END  IF 

i 

FOR  I-!  TO  Length 

Nuiner  ic  *NUM ( NumbS [  I  .  1 1 ) 

OUTPUT  Sc  USING  ”*  .  A”  : NuinbST  I .  I ) 
NEXT  I 

f 

OUTPUT  Sc  USING  "#,A"  ;CHRS< 13) 
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APPENDIX  C 


MAINFRAME  PROGRAMS 


TCAL 


TC*  L 

puwposf  :  this  poccwam  plaCcS  THFRMQCCuFi_e  c*l  i  boa  r  ion  oat 
OfcCtlwFO  FPQM  HP-SH  2  u  IC^QCQMOuTft}  INTO  * 

«0«h  iiUWK*  6  FQOM4T  ff)H  US£  IN  ThF  PPOG»A*  TF  !  X 


Dl^kNSlCN  TPPOflt (200».rAMeC200l  • T NO  2  2 (2001 
1  •  1 

l  0  Ofc  *  L>  (  0  /  .  *  I  r  AMf'  I  I  ) 

IF!  I »mM l!  J  .F0.-200. o  |  GQ  TO  20 
pt  aui  o  r  .  *  i  r  PwofjF  (  i  j 
pfc  4UI  0  r  •  5  |  TNU22  I  l  I 
n i rt  i 
I  - 1  ♦  I 

go  r  o  io 

20  CuN r i Nut 
C 

OO  JO  *  -  I  • N I  TF  MS 

•  W|  TF  lo,30  I  I  ,rAMH{  I  I  .rpOUHEl  I  I  .  r NO 22  1  t  | 
jo  con  r i nuc 

40  FUMMA  I  t  rt*  .  •  AMfll  tN  T  4C  I  •  .  2*  •  'UMfltiE  < C  »  •  . 2* .  'NO/ 2lE  fCIM 
bO  F  \'M  M4  |  |  U,|  l.ux.2tFo*J.o<)> 

STOP 
F  NO 


puwpost:  this  program  performs  a  first  u»oew  curve  Fir  «y  thp 

LEAST  saUAWfS  ME  THGO  FG«  T  mC  CLlu*nS  OF  JM*.  THF 
fiwsr  COLUMN  LISTS  VALUtS  (JF  X  ANl)  ThF  SEC.JNO  LISTS 
VALUES  OF  F(x|  W  V.  FOR  THERMOCOUPLE  CALI°H‘r  luN, 
Lf  r  A  s  f  He  A»HltNI  TF «PE w A  I  L«t  .  ANO  LET  v  =  FI  ThFR 
TMt  POOHE  GR  NUZ2LL  T  E  MP e  Q  a  T  UR  £  aS  0E5l»tD. 

DIMENSION  x|  I  SO  I  .V*  l  SO  I  .  X  SO  (  I  SOI. XV|  ISO  I.VESri  ISOI 

5U  M  x =0,0 

SUM  x  So  =  O  «  O 

SuMxv=n .0 

Su*  v  - o  .  0 

Sumnum=o  .  0 

SUMOF N5U.0 

»fc AUI  t » i I N I  TENS 

oo  io  (si ,Nir£Ms 

WFA0I7.5IXM  I.YIII 
x  SOI  l  I s  *  I  l  I *52 

AV(  I  I  =  X  |  M  SV I  l  I 
SUMX  sSOMX  Mill 

SUM X  SO  =  So* X  so*  XSG  I  l  ) 

Sum  x  v  r  *,ijm  x  vmv|  |  j 
jUMV r SU*V ♦  V  |  M 
10  CONftNUt 

*  I  rt^:FLl)»I  INI  IF  MS  | 

AS | SUMY  3 SUMXSC-SUMX OSUNX V |/ | X | T EM S« SU V *SQ- Su^ 5 ? | 

O-  f  xl  r*-_MS5SuMxY-SUMX!5SUMV  |/|  X  ITEMS-SSUMXSO-SuMXt:.?  I 
VO A N = bUM Y / X | r t * s  • 

DU  JO  I  s  1  ,  N  I  T  (-:  Mb 

VF  s  f  I  I  I  =  ti*  X  |  (|  ♦  A 

•MIIFIN.4I VOST (It.  Y ( M 

SUMNUM  -  Sumnum  *|TESTI|I-VH»WI«2 
SUMO  I-  NsSumuEN*  I  VI  I  |  -  V  ^  AP 

jo  coo r i Nut 

W‘,D:‘.U  '*  NUM  /  SUM^fc  *1 

•  Hi  r»-  |  ,  f.:  J  A  «•<  VJ 


nnnnnnnnnn 


JETCURV 


JE TCUWV 


pupposf  :  this  ppoo.ham  fit’I  jet  tr*  jfcrc«v  data  to  the  eojation 
/=ASV/|tJ*Y|  py  IMfc  least  SDOAWhS  **ElHQn  AND  c  QMF*U  r  F  a  A 
CUPHMLATfGN  CCfeFF  ICIEST  DATA  POINT*  APl  w  E  a  o 

IN  F  Ob  F  FJPMAT  h«G«  UNIT  "  7 w  C  I  SK  . 


I MfcNS ION  Y(SO»./(SO).VV|SOI./Z{SO)  .YSQISOl.YZfSOI  *ZFST ( SO  1 
So«y  =  u . o 
<iU«VS.jsO.O 

SLMV /  =  0  •  O 
5UM^:0. 0 

SUNNUMs  0 . 0 

SumoEN  =  U  •  0 

QfcAOl  f • » I N | TFMS.FACTCW 

oq  io  i = i ,n i re **s 

kE  AO  I  7 . 5  1  r  I  |  »  .  z  I  I  I 
Y  f  I  I  -  Y  |  I  l/F  »CTQP 
YY (  I  |:|  .  0/Y  f  I  I 

I  I  =  i  .o/i z  i  r  >/f actor  i 

YSu (  I  I  S  VY  f  (  ) 

YmisYYlU^iMU 

jUMYlSUMY*  T Y  |  [  | 

>uMYsa  =  sumvsu * ysu (  i  i 

SU*AY  l  =SUMY  l  *Y  Z  I  I  ) 

SUM/ :$u »z*  [  I  » 

l 0  CUN  T  l NUL 

x  1  HnI  f  F  M  s  I 

AA:<Su^/-5JUMy‘J0-SLMY5SUMY/)/C*|  T  E  MS  5  S  L  *  y  SO- S  U«  Y  S  *  4?  1 
UH  =  I  x  i  rt'^ibUMV^-sUMY^So^nM  XI  fFMS;su4'rW-SvjMYt*2  ) 

Ail . 0/ AA 

t’sAyHH 

ITEMS 

DU  4»0  I  S  I  .  N  1  T£MS 

ZFb  r  (  I  I  s  a(K:Y  Y  (  |  I  ♦  A  A 
-P|  Tclo,  JOI/ZI  I  I  • /E  ST  I  U 

SUMNUM  =  SUMNUM  ♦  I  /  t  «  r  I  l  \  -  /tiAO  |  SS2 
SJJMOL  Ni  GUMUF  N*|//|II-/bAW|$^^ 

/o  cun  r i Nut 

0b«)  =  Su^NUM/SUMOEN 

«o  t  rr  i  *,,40  1  A.u.kso 

JO  EU  *MA  f  t  I  <  ,  •  /  (ACluALl  =  •  ♦  F  I  2  .  fl  .  J*  ,  *  /  IE  jH  s  •  •  F  |  2  •  P  I 

4  0  FORM  A  |  I  l  X  .  •  A=  *.Fl?.0.2x.'^i  *  .M  ^  ,  d  ,  i  *  .F  I  /  .  tt  I 

S  l  Dm 


4.  GRAB 


c 

C  Gft*0 

C 

c  purpose:  data  transfer  fro**  the  hp-re2  fcM  i  croccmpu  ro« 

c  to  iHt  IBM. 

c 

C  HV  (JILL  CUL  tJPfc  T  m 

c  FOR  FALL.  QUARTER.  1R82 

C 

C  F  I L 1 9 b F  0*5  TFhMINal 

C  F  l L t O E F  0e>  TFRMtNAL 

C  f  ILEOh^  0/  DISK  mvdaTa  DATA  I  PERM  | 

c 

C  GLOBAL  TKTLIH  F0RTM002  MOC2EE*- 

C 

C  TYPE  IN  rwt  ABC  VE  <♦  LINES  TO  MAKE  THIS 

c  FQRfRAN  PROGRAM  RUN. 

C 

C 

dimension  OATAl 30001 
I  =  I 

AW  I  TF  <  6  .  flo  I 

80  FCJRMA  M  2X  •  •  OEG  r  N  INPUT  INC#  DATA  FROM  THE  | 

c 

c 

10  CONTINUE 

PE  AUI s . s I  OA  T  A  (  I  | 

1  =  1*1 

IFIUATAI  I-|  I  .NE.-200  I  Gore  10 

c 

N  |  TtMS  =  I  -  I 

8  FORM  A r ( 2*  I  | 5 .  •  OATA  POINTS  *E  RE  £Nff«FC.'  I 

-liiu^.oi  n  i  r  e  m  s 

c 

C  NO*  T  n  A  T  ALL  data  HAS  E  N  fc  N  TE  WF  Q .  «fi  l  Tp  IT  JUT  ON 

C  OISA. 

c 

S  FORMAT  I  2  #  •  '0*1*1  *  *15**  I  =  ••  IF  is* SI 

I  *  I 

20  *W  I  r  F  I  7 ,9  |  OA  T  A  I  I  | 

1  =  1*1 

(F  I  JA  T  A  |  f |  NC  . -20C  )  GO  T  C  20 
C 

C  ALL  data  has  BEEN  .UiriFN  ONTO  OISk. 

c 

STOP 
E  NO 
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TDATA 


roAr« 


PljRPUSE  J  Tm(S  PQOGHAM  PLACES  The  “UOVAKT  jet  data  ofcfi  VEO 
fwom  rnt  m  i  c«uc  o*Ru \f  *  into  a  «ijoe  UROFfttr 

FlJWMAI,  it  Ai.\>C  CQNVtHTS  UKITS  QF  LENGTH  FROM  bl 

TO  -FIWIC. 

01 MfcNb  CUN  X l ?00  I  , V  C  2  00  I  .  Z <2  00  \  , rPOOttE  \  2  00 )  .  T  AMB 12001  • TNO  Z  2  C  2  00 \ 

I STOt V ( 200  I  |AMM|  2001  ,tMM|  2001  • 2NM  j  2001 

I  -  t 

10  «e ao(ot,«i  *  c  i  ) 

in  *(  1  1  ,eu.-;oo.oi  GO  to  20 
**MMI=2C>. '•£*(!» 

Mi  |  U  t  .  ?  )  *  C  1  ) 

YMM |  |  |=2S.4$V|  I  » 

«t  A»)C  0  r  I  L  I  I  1 

Z-* *  I  t  >  1  > 

Rt  *Ul  0  f . i J  TP-nee  1  f  1 
Wf  AIJCO'.SI  f»MH|  1  | 
at  MM  0  f  .  *  )  1  NU22  C  I  ) 

Rfc  a  ij  (  j  7  •  i  |  S  T  UF  v  (  1  I 

N  1  T  h  M  b  :  | 

I  -  t  ♦  I 

r,u  » i)  1  u 

20  cum  r 1  Nut 

•  m  1  rr  (  s  .  ao  ) 

OU  JO  I  s  I  .  N  1  TfMS 

•  Rl  ffclt,.t)OI  f  ,(VM|  1  1  .ymm|  I  I  ,2MM(  t  »  •  T  PROSE  l  t  l  %  T  A*<M  I  I.CnO^IH 

1 s  rot  veil 

jo  cun  r  r  Nut 

40  FUUM  A  (  f  tJ*  #  •  X  (MM|*.4tf,'Y  \  M  ¥  »  •  .  2  X  .  •  P  OH  1C  I  *  .  I  *. 

i  ‘amuilnc  i c i •  .  i * • • nc /ztc  (Cl •  .t  x. 'srn  otv  iC)*t 

SO  clJRNAr(lX,|J,J*,2(P/*J#J,<i*»-7.J.4A,^CF<S.J,AJ<|| 


r»nnnnnr»  n  no  rv  n  n  n  n  n  onnnonnnrinrunnnnnnnnnnoononn 


C0NT0UR4 


con  ruuH a 


purpose:  this  0«oc154M  is  designed  rc  cispl^v  buoyant  jft 

tE^Ot«»TO«t  i)Afl  >Sl*G  A  CCNIOOP  PLOTTING  PACKAGE 
ava  i  l  a  t  «  i  r  h  oisspla.  thf  fallowing  pa  w  D  A  T  A  is 
w  f  a  i)  fhc«  a  oisk  p  iLt :  omcpf  positions  in  xyz 
COORDINATES  PECAllvF  IQ  th^  noz/lE  TIP.  AMBIENT 
f  p  mPK  P  a  r  o»e  .  PROBF  TtMP^PAru^P  «NO  nG/ZLc  TFmPFRATuRE. 
fHt_  PQOf't  ANG  NOZZLE  rE“PEF4  IUPFS  APfc  T  R  A  N  SF  •  j  P  M  F .  j 
MY  rAL  |  HP  4  f  l  QN  CiJEhP  ICICNT^  ANO  NOP  *  *t_  I  Zfc  0  WW  T  4MfitfNT 
rFMPt^ATuPfc.  Th£  XYZ  COORDINATES  APE  CONVtPTED  T  Q 
xfiw  CC3CPCINATFS  RELATIVE  TC  I  He  l  N  T  F  W  SE  C  T  |  QN  QF  T  Hfc 
u a  r  4  plane  «i  in  rnt  CENTt.PL  i  ne  trajectory  qf  thf  jet 
ithf  s-4x  ib  i s  tangfnt  rn  iHt  tr  a  jtc  top  y  j .  ihi  5 
PPUGWAH  ALSO  CCMPWTFS  IHt  «4  IF  QF  hF A  r  TWANSFEP  FWUM 
IHt  Jfc T  rc  THE  AMHIENT. 

NOTE:  TpE  FOLLOWING  VAfcuFS  MUST  BE  INSERTED  as  the 

FlwST  l  I  nE  of  cata  in  FPfcE  fcpmai:  thf  total  number 
QF  OAT  A  POINTS  IN  ITEMS!.  TkF  ACUTE  4NClE  EET^EFn 
IHt  Q4  14  PLANP  ANO  HORIZONTAL  (  Thf  1  A  |  AM)  THF  VERTICAL 
DISTANCE  °E  T  a  c  t  N  The  Z-AXIS  ANO  the  INTERSECTION 
POINT  ClSCbSSEC  ABOVE  |ZA) .  THfc  CENTERLINE  VELOCITY 
IN  M/s  ivEL)  AND  JfcT  WIDTH  | N  MM  (WIDTH). 


dimension  X  (  1 00  I  .  Y I  I  00  »  .  r  I  100  I  .  TP  |  100).  TA  I  100).  TN|  100  )  .  TMA  T  |  10.  10| 
COMMON  «OhK| I6D00I 


-  -  -  -  HEAD  DATA  FROM  UNIT  "7"  0  f  SK  *  $  $  5 


pfc  A  o I  7  .  $  1  NfTEMS.THErA.ZA.vEL.wICTH 
DO  20  1=1. N ITEMS 

WE  AO  I  7  .  *  I  4,X|  !  I.A.YI  ll.TPHI  .TAl  (  )  •  T  N  (  I  )  .  A 
y  {  I  l  =  I  Y  {  I  )  -z  A  |  /SIM  The  r  a  ) 

CALIBRATION  COEFFICIENTS 
TPI  l  Jr.y*aeOQp7SSTP(  I  I  ♦  I  .  7*0  7RbCS 
TNI  l  )=.  M  322  H2eiIN|  I  )  ♦  I  .dWSAOeb 

rill  =  (TRIM  -  TAlili/ITNlll  -  T  A |  I  )) 

2  0  CON  T  I NUt 

3  0  Forma ri/K,  •  I . a . y  .  r  =  •  ,  I  S  .  IF  I b . 5  ) 


i  5  5  ;  U  I  ms:  N  S  I  0  N  QF  TMAf 
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c 

c 

c  ~ 

c 

c 


c 


I xn i m  s  io 
I  V!)  |  M  ;  10 


■i  j  t-  I  NO  Thf  maximum  $nq  m|N|*um  r  fc  m*jf  h  a  TURF  3  S  S  ^ 


IMIN  =  tOO.J 

r  M  A  X  r  (J  .  J 

x  m  a  X  -  J  .  O 

<  M  1  N  =  I  O  u  •  0 

Y M  I  N  =  10  0.0 


QU  A  il  |  s|  ,NlTt«b 

I  F  I  T  |  I  |  .U  .  T  m  a  x  I  f  m  a  x  =  T  (  M 
I  F  l  M  I  )  .Li.  iv|M  pin  =  TMI 


IFIXJI  |  .of.ivun^UMM  i 

IF  I  ■  4  t  )  •  I  *  ,  l  «  |  M  I  "*  |  N  :  I  I!  i 
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CUN  I  |  M>c 
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^UAOJwANf,  T»-.P  ISt  AXIS  IS  Cf  NTc^eO 
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Cr--|  |  Nf/fLtMri  > / »  000  . 
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A  N  fe  A  =  i;  A  5:0  Y 

O  ESTABLISH  TMfc  [NC«E*fcNT  SIZE  FC3  CONTOUW  PLOTS  *  z  z  z 
r I NC«= I TMAX-TMjNI/5,25 


z  GENERATE  fHE  GWIC  z  s  3  s 


CALL  COMPHS 
CAi  L  TEK6  1  e 
CALL  pale  I  d  •  •  ^  •  I 
C  At  L  IJLU«UP  10.  /5I 
CALL  PHY  SuP (  I  •  Z  S  •  1  •  ) 

call  A«t azo \ s . s  . s . s  i 

C  all  HFIGhTIO.IdOl 
CALL  C  AW  TIJC 
Call  i n  r  a  x  s 

CALL  xna«l(»x  I«WI4*.100I 
CALL  YNA-fct**  (M»*hMOCI 

c  all  *  r i c«  s ( 2  i 
c  all  v  r  i  c«  s  i  4?  i 

CALL  Y  AX  ANU  I  ‘iO.OI 


C  all  S  *  I  SS  * 

CALL  MEAOI  NI  •  TFMOtO  A  f^O(=  CONTOURS  IN  A  tuijYAN  T  jF  T  1  *  ,  100,2. 

C**LL  mFAOINI  •  *|  th  a  r  R  0  S  Sf  l  U  «•  l  N  G  *Mri  I  EN  '  •  (  00  •  2  •  •  4  | 

CALL  hFAO|N(  •  (  2  S  X  PLC* 

CALL  « rl  a  o  |  N  |  •  l  F»  L  A  N  E  A  |  %  •  ,  |  o  0  •  I  .  .  A  I 

CALL  GO AF  ( -  30 • 0 . OS . d  .  JO . 0  •-  JO .O.Ob.O.  20.01 

-  GENERATE  ThF  I  N  r  e  PPCL  a  f  E  U  TFMOtRATLRc  *  a  T  p  |  x  "f^Af  *  5 

CALL  clCi)*(JN  I  1  oOOO  | 

CALL  ZP ASL  (  f -  I N  | 

CALL  f»GN*  A  f  I  1  X  0  I  M  ,  I  VO  I  *  I 

CALL  Gt  Tm*  M  X,  Y  ,  I  I  IF  ws,0  I 
call  FM)«*r  i  r«Ar  ,i)i 

-  iOLVt  THt  PATE  OF  MEAT  TW»NSFen  A  5  9  S 

jSO«=0..1 

GO  r 'j  I  =  |  •  l  x  u  i  m 

nu  •> «o  j:  i  ,  i  »ni« 

R»1  =  (  I  .  o  ■)  U  *)*  !  J  J  .  i  I  P  I  >F  -C,£  r  *  A  I  (  1  .  J  | 

-*>.  {  !  ,  J  |  ~  C 

*  J.<.U''lr,<>}F.»c.rvA|  (  i  ,  j|;  *  j 

-  GNOr>/f.--|  |il  U!  4  I  »  j  |$$«  |»  lOu'). 

M:iJU«U<-FliaI  I  )-[  |  *  fj  K 

Y  V  SIJvjAl.V  -FUlAri  I  -  I  J  :  L  Y 
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APPENDIX  D 


TABULATED  DATA 
TABLE  1 


ROTOMETER  CALIBRATION 

ml /s  Std  Dev  ( m  1 , '  s  ) 


5.94 

0.32 

8.02 

0  .  56 

10.05 

0.41 

11.85 

0 . 24 

14 . 03 

0  .  33 

15.98 

0  .  58 

17.41 

0 . 31 

19.09 

0.34 

20.75 

0 .23 

22 . 38 

0.28 

24 . 27 

0.22 

26.03 

0.20 

27 . 64 

0.37 

29 . 53 


0.24 


TABLE  2 


TEST  RESULTS 


Crossflow  velocity: 

Nozzle  flow  rate: 

Nozzle  inside  diameter: 

Nozzle  discharge  velocity  (mean) 
Nozzle  temperature  (mean) : 
Ambient  temperature  (mean) : 
Froude  number; 

Michaelis-Menter  equation: 

Coefficient  A: 
Coefficient  B: 


.13  m/s 

11.85  ml/s  (25% 
7.144  mm 
29.558  mm/s 

41.8  °C 

24.9  °C 
14 . 8 

2.64284325 

1.03698254 


Plane 


Y  (mm) 


4>  (degrees) 


7.327 

46 

1.20 

21. 370 

58 

6.392 

39.688 

70 

15.349 

62.889 

80 

18.^83 

87.313 

86 

27 . 62d 

[•] 
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